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A STIRLING RNGINE COMPUTER MODEL FOR 


PERFORMANCE CALCULATIONS 
by Roy Tew, Kent Jefferies, and David Miao 
Lewis Research Center 

SUMMARY 

To sumx>rt the development of the Stirling engine as a possible alter- 
native to the automobile spark -ignition engine, the thermodynamic char- 
acteristics of the Stirling engine were analyzed and modeled on a computer. 
The modeling techniques used are presented. The performance of an exist- 
ing rhombic -drive Stirling engine was simulated by use of this computer 
program, and some typical results are presented. Engine tests are planned 
in order to evaluate this theoretical model. 


INTRODUCTION 

The Department of Energy (DOE) has established prc^rams whose pur- 
pose is to reduce fuel consumption and emissions of highway vehicles. The 
Stirling Engine Highway Vehicle Systems Program is one such program. 

Its purpose is to develop the Stirling engine as a possible alternative to the 
spark -ignition engine. It will be implemented through government participa- 
tion with industry. NASA Lewis Research Center has project management 
responsibility for the program. 

At NASA Lewis a Stirling -engine digital -computer model is being devel- 
oped for predicting engine performance. This report documents the present 
modeling techniques and shows some preliminary predictions of one version 
of the mcdel, which is configured to represent a particular single -cylinder 
rhombic -drive engine. 
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Analysis of the ideal Stirling cycle is straightforward (ref. 1), but this 
analysis is useful only as a reforence against which to compare real cycles. 
Sdunidt derelqped a somewhat more realistic analysis that considers the 
^Ict of dead volume and assumes sinusoidal piston m(Aion {rei. 1). The 
basic power and efficiency ^e Carnot efficiency) calculated by the Schmidt 
analsrsis can be multiplied by experience factors to estimate performance 
in tlw initial stages of engine sixing. (Ranges of eiqierience factors are 
discussed in refs. 2 and 3.) Rios (ref. 4), Qvale (ref. 5), and Martini 
(ref. 2) have eadi developed models that first calculate a basic power and 
efficiency and then correct these basic values with sq[>arate, indqiiendent 
loss calcolations. The most detailed performance models so far reported 
simulate real-time variations in gas temperatures, pressures, and flow 
rates at various control volumes within the working space. Some or all of 
the loss calculations are an integral part of the model and can thus internet 
with the basic thermodynamic calculations and with each other. Detailed 
models that have been r^x>rted in the <^n literature are those of Urleli 
(refs. 6 and 7), Fhriielstein (ref. 8), and Finegold and Vanderbrug (r^. 3). 
A listing of Uriel! *s computer model is given in reference 8 and a listing of 
the Finegold -Vanderbrug model, in r^erence 3. Allan Shock of Fairchild 
Industries is also developing a detailed model (ref. 9). . The more detailed 
modtis should be oi value in studying the effects on performance of ei^ne 
details not considered in the simpler models, in refining the loss calcula- 
tions of the simpler models, and in investigating the significance of phys- 
ical effects such as gas inertia and pressure dynamics on engine perform- 
ance. However, these models appear to require too much computer time 
for use in a design optimization program: One run with .such an optimiza- 
tion program might investigate hundreds of designs, or even several 
thousand. Since the less detailed models such as those of Qvale, Rios, 
and Martini require much less computer time, they would seem to be more 
^qpropriate for use in a design optimization pre^ram. Unforhinately, the 
models that have been published in the open literature have not been ade- 
quately validated. 
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N. V. Philips snd Us llconss«s have both a design -optimization pro- 
gram and detailed performance prediction models* which have been devel- 
oped over the years with the aid of extensive test data. However* these 
programs and the test data are proprietary and are presently available 
only under license. 

The model being developed by NASA Lewis differs from the Uriell 
and Finegold models essentially in that the Lewis model assume.s that gas 
inertia and pressure -wave dynamics can be neglected in predicting Stirling 
engine performance. The Lewis model also uses an integration technique 
that avoids instidoility caused by excessive heat tran.sfer between gas and 
metal when large time liu:rements are used. Thus, the Lewis type model 
should be inherently more efficient in terms of computing time, than the 
more general Uriel! and Finegold models. 

The Lewis model differs from the Hlos and Martini nuxiels essentially 
in that it more closely repre.sents the distributed -parameter nature of the 
working space by dividing each heat exclumger Into several control volumes 
and, also, by making the heat -exchanger inefficiencies an integral part of 
the cycle calculations. The Lewis nuxtel is thus more general in nature but 
less efficient in terms of computing time than the nuxlels of Kios , .id 
Martini. 

Two versions of the computer performance nuxtel have been devoloptxi. 
One version is configured to model a rhombic -drive grouixl pi>wer unit 
(GPU) designed and Ixillt for the U.S. Army by General Motors. This unit 
was designed to prtxluce 3 kilowatts of electric pt>wer (or about t> kW of 
br;ikc piiwer). The other version is conflgureti to nuxlel a free-pl.ston en- 
gine. This doi'uinent contains a description aixl program listing of tiie 
GPU version of tlie model. The GPU engine i>arameters are given in 
refeivnce 2. The model is also briefly de8crllH>d in reference 2. along 
with some performance prediction comparisons between it and the Miu'iini 
nuxtel. 

Tlie GPU -3 engine has been run thrvmgh initial checkout tests at N.ASA 
l.owis; engine accessories were powortxt by the engine and instrumentation 
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was minimal (ref. 10). Testing with accessories run independently of the 
engine and with more instrumeirtation is now getting under way; these tMts 
will provide basic engine data for evaluation and refinement of the computer 
model. 


GENERAL DESCRIPTION OF MODE L 

The GPU engine, for which this model predicts performance, is shown 
schematically in figure 1. There are actually rS^ht separate flowpaths 
through each of eight regenerators and coolers with five heater tubes serv- 
ing each regenerator. In the model, however, it is assumed that the same 
flow conditions exist in each of the eight flowpaths so that it is necessary to 
model only one path. The model represents the working space by a series 
of subdivisions called contrc'. volumes; this type of model is sometimes 
called a nodal model. 

The model calculates indicated pv>wer and efficiency for given engine 
speed, mean pressure, and fixed heater and cooler metal temperatures. 

The indicated efficiency is based on heat into the gas plus conduction losses. 
The model also simulates temperature, pressure, .and flow variations over 
the cycle at vai'iuus stations in the working space. The working space con- 
sists of the expansion space, the heater, the regenerator, the cooler, and the 
compression space. 

The engine working space is represented by 13 control volumes, as 
shown in figure 2; the adjacent metal walls are represented by 13 corre- 
sponding control volumes. The metal temperatures, except for those in 
the regenerator, are assumed to be constant. This is a reasonable assump- 
tion for any given run since the heater and coolci metal temperatures are 
essentially boundai'y temperatures that are controlled by the combustor and 
the cooling water flow, respectively; these temperatures vary little over a 
cycle because the metal heat capacity is much greater th.an that of the gas. 
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Th« calculation procedure used in the model Is outlined iu figure 3; 
the equations used in the calculations are discussed in the section EQUA- 
TIQNS AND ASSUMPTIONS US E D IN DEVELOPING THE M ODEL . Each 
set of calculations shown in figure 3 within the inner loop is made at each 
integration time step during each engine cycle except that it is necessary 
to use the pressure-drop equations only during the last cycle. Also, the 
conduction and shuttle heat -transfer losses are calculated Just once, dur- 
ing the last engine cycle. (Shuttle heat transfer Is heat transfer by heat - 
ing of the displacer at the hot end of the stroke and cooling of the displacer 
at the cold end (ref. 2). ) Between cycles, corrections to the regenerator 
metal temperatures ai*e made to speed up convergence to steady operation. 
Typically it takes id>out 10 cycles with regenerator temperatui'c correction 
to adjust the regeneridor metal temperatures to their steady operating 
values. (The model predicts that these temperatures vary with an ampli- 
tude of about 6® C or less over a cycle. ) In addition, a number of cycles 
are required for the leakage between the working aiKl buffer .‘’.paces to .ad - 
Just the mass distribution. The smaller the leakage rate, the longer the 
time required for the mass distribution to reach steady state. For the 
range of leakage rates considered thus f;u\ it takes longer for the mass 
distribution to steady out than for th.^ regenerator metal temperatures tv) 
settle out. Current procedure is to turn the metal temperature conver- 
gence scheme on at the fifth cycle and viff at the 15th cycle. The model is 
then allowed to run for 15 to 25 more cycles to allow the mass distrlbutlv)n 
to settle out. When a sufficient number of cycles have been completed for 
steady -state operation to be achievevi, the run is terminated. 

Current computing time is about 5 minutes for 50 cycles on a UNIVAC 

1100, or 0. 1 minute per cycle. This is based on 1000 iterations per cvcle, 

-5 

or a time Increment of 2M0 second when the engine frequency is 50 hertz. 

The computing time can be decreased by decreasing the number v>f U - 
orations per cycle; the corresponding effect on predicted indicated power 
and efficiency is shown in figures 4 uikI 5 for i>ne series v)f runs. Tliese 
figures show that .some error is iiu oduced by reducing from 1000 t») 500 
iterations per cycle and that the error becomes more significant if v)nly 
200 iterations per cycle are used. 
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Th* coirr'rter prograni is written In FORTRAN V and, In card format, 
is about 1300 cards long (including plotting subroutines). The program was 
designed to be an engineering tool for use in est^lishing the validity of the 
modeling techniques. It could be used rather easily by others but has not 
been designed or extensively documented for that purpose. A listing, not 
including the plotting subroutines, is presented in sq>pendix E. 


EQUATIONS AND ASSUMPTICtt^S USED 
IN DEVELOPING THE MODEL 

' ' ' y ' ■ 

First, the basic equation.^ and assumptions used in nuUditg the thermo* 
dynamic calculations in the working space ai'e stated. Then a relatively 
coniplete presentation of the equations used in the nuKiel is made that cor- 
responds very closely to the steps simwn in the outline of calculation pro- 
cedure (fig. 3). 


WORKING -SPACE THERMODYNAMIC CALCULATIONS 

Each of the 13 gas control volumes shown in figure 2 is a special case 
of the generalised control volume shown in figure t>. The generalised con- 
trol volume includes flow acrv>ss twi> surfaces, heat transfer across a sur- 
face, and work interchange between the gas and a piston. Each of the three 
heater, five regenerator, and three cooler control volumes has flow across 
two surfaces and heat transfer across one surface Init is of fixed volume; 
therefore, no work is done by the gas in these volumes. The expansion - 
and compression -space control volunies each have flow across one surface 
and heat traiisfer across one surface and ;u'e of variable volume; there- 
fore, the gas in these twi» volunu's is respivisible for the w\>rk output of 
the engine. 
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The three basic equations used to model the thermodynamics of the ga;? 
in each control volume are conservation of energy, conservation of mass, 
and equation of state. These equations are used to determine the temper - 
atiure and mass distributions and the pressxire level within the working space 
at a particular time. A fourth basic equation, the momentum equation, in 
steady -state form, is used to calculate pressure drop across each control 
volume in order to evaluate its effect on indicated power and efficiency. 
However, this pressure -drqp calculation is decoupled from the thermody- 
namic calculations; it has no effect on the temperature and mass distribu- 
tions. This assumption simplifies the model and should be reasonable when 
the pressure drop is sufiiciently small in relation to the pressure level. 

The energy, mass, and state equations, as written for the generalized 
control volume shown in figiure 6, are as follows, where three formulations 
of the equation of state are shown: 

Conservation of energy (for negligible change in kinetic energy across the 
control volume); 


hA(T^ - T) + Cp(w^Ti - w^T^) - 


dt 


( 1 ) 


— (MC^T) 

Rate of change 
of internal energy 
of control volume 


Rate of 
heat transfer 
across boundary 
of control volume 


Rate of enthalpv 
flov across 
boundai y of 
control volume 


Rate of work 
done by gas in 
control volume 


Conservation of mass: 


dt 


= Wi - Wo 


( 2 ) 
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Sqa«Uon of 


PV » MBT for IdMi g«s 


PV « MR |t ■»> M558 pj for hydrog« . - real gms 

PV * MR |t ^ ^.OieiS for heUttin • real gas 






where 

A heat transfer area of coidrol volwne 

Cp,Cy heat capacities at conr^ant pressure and volume 
h heat-trmisfer coefficient 

M mass of gas In volume 

P pressure 

R gas constant 

T bulk or aver<)ge temperature of gas In volun^e 

temperatures of gas (lowing across surfaces I and o, 

respectively (In tig. 6) 

temperature of metal wail adjacent to heat -transfer .o'ea A 
t time 

V volume 

W|,w^ flow rate across surfaces I and o, respectively 

(The real -gas equations of state were developed from data In ref. 11. ) 
Several assumptions are inherent In the use of these equations: 

(1) Flow is one dimensUmal. 

(3) Heat conductUm through the gas aiKl the regeneratiu' matrix along 
the flow axis Is neglected. The thermal coiKluctlvUy of the regenerator 
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matrix is assumed to be Infinite in calculatiitg tho overaU gas-to-matrlx 
heat -transfer coefficient. 

(3) Kinetic energy can be neglected in the energy equation. 

(4) The pressure -drop calculation (based on the conservation -of - 
momentum equation) can be decoupled from the <^er three basic equa- 
tions. This implies use of a uniform pressure level throughout tho work- 
ing space at a given time in at^lying these equations. 

(5) The time derivative term in the momentum equation is neglected 
(see appendix D). 

(The last three assumptions are not made in the generalized models of 
UrleU (refs, 6 and 7) and Shook (ref. 9).) Thus diesc two models pi*ovide 
a me.ans of checking the validity of these assmnptlcms. ) 

Jr appendix A it is shown that equations (1) and (2) aixl the ideal -gas 
equation of state can be used to derive the following differential equation: 


MC ^ - hACT - T) C Wj(Ti - T) - C -T)^V^ (4) 

dt ^ dt 

The same result is obtained if either t»f the real -gas equations of state 
are used in the derivation. This equation s.iy.>^ tltat the bulk or average 
gas temperature of a control volume is a fuitotion of the following three 
pt\K esses; 

(1) Heat tr;msfer across the boundiu-y from the wall 

(2) Gas flow across the boundiu'v 

(3) Pres.surc level 

One approach to numerically Integrating equation (4) Is to solve for 
the temperature derivative 


dt MCp 


(T^ - T) + J (Tj ~ T) - -1' iT^, V.T) 
M M 


V dP 
MC dt 

r' 


(5) 
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and then set 


jUM = T* ^ — At 
dt 


where 

new gas temperature (at time t At) 
old gas temperature (at time t) 

At time increment 

This is a first-order numerical integration. A similai* but higher order 
technique (Runge-Kutta fourth orde* , for example) could also be used on 
equation (5). This type of approach was not used. 

The approach used (which helps to avoid numerical instability prol)- 
lems) was to decouple the three processes that contribute to the temper- 
ature change and solve for the temperature change due to each process 
separately. This second approach *vTows a trade-off between cotnputii^ 
time and accuracy of solution (as indicated in figs. 4 aixl 5) with much 
less concern for numerical ini.tabllities. The approach is suggesttxi by 
representation of equation (5) in the followiiig form; 


= 


dT 
+ — 




total 

due to 

dt j 

due to 

dt , 

due to cl\;u\ge 


heat transfer mixing in pressxu*e 


where 


dt 


V dP 
MC dt 


(5a) 


due to change 
in pressure 
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WjfTj - T) - Wq(Tq - T) 

due to ** 

mixing 

dt 

heat transfer 


due to 


hA 

MC„ 


fr„ -T) 



(5b) 


(5c) 


fiat ^qpp^idix B it is shown that equations (5a) and (Sc) can be integrated 
in closed form and tiiat equation (Sb) can be numerically integrated. Wh«i 
the results of appendix B are modified slightly to show just how they are 
used In the model, the resulting e^[>re8sions are 


mt-fAt 



'pt+AtV’'-"''’' 
< P* / 


ioa’) 


rpt-l-At 

*PM 




M 


t+At 


(5b’) 


rpUAt q,t+At . 

* “ *PM + 


K- 


pt+At 

PM 


^t+At^/j^+AtCp) 

) l.e 


(5c') 


where the superscripts t and t-i-At denote values of the variables at 
times t and t -i- At. The subscript P denotes the value of the temper - 
ature after it has been updated for the effect of change in pressure. The 
subscript PM denotes the value of the temperature after it ha.^* been up- 
dated for the effects of change in pressure and mixing. No subscript (as 
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on ttw l«fl tldo of oovinHon donotos th<» v«au« of tho t•mporature odor 
it hill boon updntfd for nil thrto eftoctn • chnn^o in pronnurt, mlxini» and 
float Iranaftr to or from tho matai. 

PRK8KNTATION OF RQUAT10N8 IN ORDER 
or CALCULATION PROCEDURE 

Tho oqiiaticmn considarad ao tar hava baan darivad and diacoaaad with 
rtforanca to IKa ftanaratiiad control votuma of Ripira 6. In the coniputar 
modid thasa aquationa ara appliad to aach of tha 19 control volumaii ahown 
in (ignra I, Than tamparatnraa« nianaaa* haat-lranafar eoafticiaida, flow 
rataa^ ale. « ara all aidiacriptad with an Index. The index variea from 1 
to 19 for variablea that are avara«:aa of the control vtUumea and from 1 to 
19 hir valuaa at tha iidarfacaa i>atwean control vtdomaa. The nnmimring 
procadiura viaad tor control vv>lume and interface indexes in defined in tip- 
nra 9. The equations diacueaed in thia aection include theae iiHlexea, 

The preaentattmi of the equationa foilowa the atepa ahown in the outline of 
calculation procedure in thpire 9. 

Preaaure (Btm> 3 in fig, 9). - The in'eaaure P la calculated by 



Where 

P preaaure 
R paa conatant 



IS 


li^ mwu In voinmt 
Tj nv«rai« temperakorn in ^ volnmn 
Vj ^ volume 

I an index denoting which control volume Is under consideration 
Equation (6) is obtained by summing the ideal -gas equation 

PVj«MjRTj I » 1,13 

over each of ttie IS control volumes (remembering that pressure is as- 
sumed to be ttie same in all control volumes) and then solving for P. 

If the real 'gas equation of state for hydrogen is used then and the same 
procedure is followed, the result is 


pt+At 


ti 


15 

' ' vJ - 0. 02358 

W 


1=1 


M 


(7) 


Equaticms (6) and (7) are both included in the program. Also an equivalent 
real -gas equation for helium is included. An index in the input data speci- 
fies whether a real or ideal equation Is to be used. 

Update temperatures for effect of change in pressure (step 3 in 
fig. 3) . • It was not necessary to use the subscripts P ;md PM of equa- 
tions (Sa'),(5b'). and (Sc*) in the computer nuxlel. Applying the three 
equations In sequ'^nce without the subscripts produced the same effect us 
if the subscripts had been used. Therefore, if the subscript P Is droppoit 
from equation (5a*) and the control volume index I is InttHHluced, the re- 
sult is 



H 


ipt 

M • M 
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This equation is comm<»dy used to relate t<onperature and pressure for an 
adiabatic fixed -mass process. 

Mass distribution (step 4 in fig, 8) . - The equations for mass distribu- 
tion are derived by assuroliq; that the mass redistributes itself in accord- 
ance with the new volun\es and temperatures in such a way that pressure is 
uniform througthcnd the woricing space. This assumption, of course, in- 
troduces inaccuracies in the solution, but the inaccuracies are rtiatively 
small because the pressure drop in Stirling engines is usually a small por- 
tion of the total pressure of the fluid. The pressure P throughout the 
working space is derived from the perfect -gas law as follows: The perfect- 
gas law for the control volume can be written 



PV, 


RT, 


Summing over the i3 control volumes 



and solving for P/R gives 


^otal 
13 V 


P 

R 



IS 


Now sidMititiilliig lOr P/R into the perfect-gas e<iiiiMon for the con 
trol vuiunie gives 

M,.ilS25aL!i 

The form of this eqastimi used in the model is 



(where represents Tj updated for pressure but not for mixing and 

heat tranmer). 

Th* preceding equation calculates the new mass distribution for the 
case 01 a perfect gas. The following equation, which can be derived In the 
same manner, Is used to ai^proximate the real properties of hydrogen: 




le 

A similar equation that aMuroximatta the real propartias of halium is tn> 
eluded In the modal. 

Flow rates (step 8 in fig. 3) . - Once tlia new mass distribution is 
known, the new now rates are ealeulatad from the old and new mass dis- 
tributions aceordinp to 






At 




and 





(or I • a. la 


o 


m) 


where w^ is the flow rata at the intartaca between eontrol volumes. 

Update temperatures in each control volume %r effect of gas flow 
between control volumes (step 6 In fig. 3) . • The followii^ equations 
(modifleations of equation (5b*) in the section WORKIKQ -SPACE THEKMO- 
DYNAMIC CALCULATIONS) were used to update temperature for the mix - 
i»ig effect followiiqc gas flow between control volumes: 


.t^At 

l.PM 




M 


t>At 


,^t..t^t /^t^t.t^t , *t+At l^tV. 

UAt .“l‘l.pn^-l "*^l , . . .. 

I pj^ * « a 1 — for 1 « a. 12 

^.At 

. ^is^is.p^ria ^12 


(la) 


M 


tfAt 

13 


mt'^At 
*13. PM 
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(where the subscript PM indicates the temperature has been updated for 
pressure change and mixing). 

The equations for the first and thirteenth control volumes (expansion 
and compression spaces, respectively) are simpler In form than those for 
the other control volumes because there Is flow across only one surface 
In each of these volumes. (The leakage flow between compression and 
buffer spaces la handled Independently; It does not appear In eq. (12). ) 

The temperature of the fluid flowing across the Interface has been given 
a new variable name $ to better distinguish It from the average control 
volume temperature T and to keep the subscripts as simple as possible. 
The procedure used to update the temperature e for each Interface Is 
now defined. 

The temperature of the fluid flowing across a control volume boundary 
Is lust the bulk temperature of the control volume from which the fluid 
came > for flow from the expansion -space, heater, cooler, or 
compression -space control volumes. This is a reasonable assumption 
for these volumes since the actual temperature gradient across each is 
expected to be relatively small. In the five -control -volume regenerator, 
however, the temperature gradient Is not small. One option would be to 
Increase the number of conti'ol volumes In the regenerator. However, to 
save computing time, an alternative aj^oach was used. It was assumed 
that a temperature gradient existed across each volume In the regenerator. 
The magnitude of the gradient was assumed to be equal to the correspond- 
ing regenerator metal gradient. 

A schematic of a regenerator control volume Is shown in figure 7(a). 
Flow across both Intv rfaces Is, for now, assumed to be in the direction 
shown (which is defined to be the positive flow direction). The cross- 
hatched area represents the portion of the fluid that will flow across in- 
terface I during the time increment At. The assumed temperature pro- 
file of the control volume is characterized in figure 7(b). The vertical 
dashed line in figxu'e 7(b) defines the temperature at the loft boundary of 
the fluid that will flow across interface I during At, If Tj is defined 
as the average temperature of control volume 1 and ATj equals 
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ont*half the changt in ttmperature across tht control volume, then 
Tj - ATj Is the ttmpersfturs of the fluid at inisrfsoe I and 




W| M 

UT 


2 ATj 


is the temperature of fluid at the vertical dashed line. (Figure % shows 
the numbering schemes used for the control volumes and the interfaces 
between control volumes, ) 

Now the temperature oi the fluid that flows across an Interface dur- 
ing At is assumed to be equal to the average temperature of that fluid 
before It crosses the Interface. The average temperature of the fluid 
in the crosshidched area of figure 7(a) is then 


1 

2 


(t, - at, . a AT,) . (t, . AT, )] - T, - AT, . lli* AT, 


Therefore, for the flow directions shown in figure 7(a), the updated 
temperatures of the fluid that crosses the interfaces during At are 


At 


^-►At ^ jt+At . ^ ^ Wj^' > 0 

* »#t 


M 


«,t+At .A 

,,t'*‘At _ rpt-fAt aTi l*l ATT* 

" M-l,P ^ 


M 


t 

I-l 


M n-l ^ 


( 13 ) 


If the flow direction is reversed at both interfaces, then 
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^ . Arr I A»r ^t+At ^ fi 

^ “M+i + <® 


^+i 


(14) 


^t^At aA 

^I-f^ = + ATj 4- ATj ^tf ^ < 0 

m{ 


Heat "transfer coefficients (step 7 In fig, 8) . - The assumptions and 
equations used in calculating heat-^ansfer coefficients are discussed in 
appendix C. The heater and cooler equidlons are based on well- 
established) steady -flow correlations for tubes. The regenerator equa- 
tion is based on an extrtqx)lation of a steady -flow correlation. There is 
a need for additional steady -flow heat -transfer data for Stirling engine 
heat -exchanger components, especially regenerators. In addition, data 
are needed to determine how to modify steady -flow correlations for the 
periodic -flow conditions that exibi In Stirling engines. 

Update temperature in each gas control volume for effect of heat 
transfer between gas and metal (and determine heat transfer between 
gas and metal) (step 8 in fig. 3) . - This temperature vqxlate is accom - 
pushed by using the following equation (a modification of eq. (5c*); 




I 


I, PM 


Lt ™,t+At \ 
\^W,I ■ M,PM/ 


1 - e 




L 


I = 1,13 


where T^ j is the wall temperature of 1^ control volume. Note that, 
no matter (low large the heat -transfer coefficient, the gas temperature 
cannot change more than the AT between the wall and the gas. Thus 
this calculation cannot cause the solution to become unstable, but it can 
lead to significant inaccuracies if the time increment At is made too 
large. 
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The heat transferred between gas <uid metal Is then calculated from 

Cp I » 1 , 13 (15) 

SO that heat transfer from gas to metal is deftned to be positive. 

Regenerator metal temperature (step 9 in fig. 3). - The equation used 
to update the metal temperatures in the five regenerator control volumes 
is 


M,C— 1 = 5,9 (16) 

‘ dt * 

where Q| is the rate of heat transfer between gas and metal. This is in- 
tegrated numerically by setting 


w,l 







MjC 


At 


(17) 


where 

M| mass of metal in volume 
C thermal ciquicitance of metal 
At time increment 

For most regenerators the thermal capacitance of the metal is so much 
larger than the thermal c^acitance of the adjacent gas volume that an me- 
cessive number of engines cycles (from the point of view of computing time) 
are required for the metal temperatures to reach steady state. Therefore, 
it is necessary to apply a correction to the metal temperatures after each 
cycle to speed up convergence. The method used is discussed in the sec- 
tion Convergence scheme for regenerator metal temperatures (step 19 in 
fig. 3). 
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Pr»«8ur»-drop calculations (st^ 10 in fig. 3) > - Since the pressure - 
drop calculallons have been decoupled from die heat- and mass -transfer 
ealoulatlonS) pressure drop needs to be calculated only over the last cycle. 
The indicated work calculation can then be corrected for pressure -drop 
loss. 

A general form of the conservidion of momentum equation for one- 
dimensional flow is 


^(pv) 

Rate of 
accumulation 
of momentum 
per unit 
vtdiume 


- JL (pv^) 


Rate of 
momentum 
gain by 
convection 
per unit volume 




Rate of 

momentum gain 
by viscous 
transport (flic- 
ttonal forces) 
per unit volume 


Ox 


(18) 


Rate of 

momentum gain 
due to pressure 
force per unit 
volume 


where 

p density 

V velocity of flow 

f friction factor 

tvvdraulic diameter 
P pressure 
t time 


X distance 

m appendix D it is shviwn that by combining the Cimtinuity iuid monu'utvuu 

equations and then noglooting tho time dorivativo torm iix tho ro- 
sulting equativvi, the Mlowliig oquatUai results: 


V di* 



dx 4 


dP 


0 




{\ 9 ) 



This equation can be integrated over a length L for the special cases of 
adiabatic or isothermal flow processes (the two extremes). When the re- 
sulting adiabatic and isothermal expressions are ^lied to the GPU re- 
generator, the contribution of the v dv term is negligible for the two ex- 
tremes. Since the effect of the term is more significant in the regenerator 
than in the heater and cooler, the expressi(xi for pressure drop can be re- 
duced to 






( 20 ) 


or iq)plying the differential equation (20) over a finite length L 


AP = i-ipv^L (21) 

where AP is the pressure drop over length L. 

A modification of this equation can also be used to ccount for the ef- 
fect of expansions and contractions in flow area. The form of the modified 
equation is 


AP=Kipv2 (22) 

2 

R is applied at each area change in the flowpath between the expansion and 
compression spaces. At a particular point where an area change occurs, 
K is a function of the two areas and the direction of flow (since an expan- 
sion for one flow direction is a contraction when the flow reverses). The 
term K is calculated in accordance with the procedure given in refer - 
ences 12 and 13. 

For the heater and cooler control volumes the friction factor f is 
deter mii ed from 
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*1 = 




16 


N 


Rej 


0.046 


N 


0.2 

ROj 


N„ < 1500 1 




(23) 




J 


where is the Reynolds number (based on plots of smooth -hibe fric- 

tion factors in ref. 14). 

For the regenerator control volumes 


-O.OlOONp 

fj = 0.96e ^ + 0.54 1 = 5,9 (24) 

which is a fit of the curve shown in figure 8. This curve was derived from 
the experimental steady -flow air data taken on the regenerator cooler unit 
of the GPU -3 engine and was extrapolated to the low Reynolds number 
range (below 50). Since the sample runs for this report were made, this 
cxirve has been compared with other wire -screen data from Kays and 
London (ref. 13). This comparison suggests that the curve should be some- 
what higher in the below-50 Reynolds number range. (This range is im - 
portant in the regenerator. ) 

There is a need for additional pressure drop data for Stirlii^ engine 
heat -exchanger components, especially for regenerators. Also additional 
data are needed to determine how to modify the steady -flow pressure drop 
correlations for periodic -flow conditions. With the pressure level P 
known (which is now assumed to be the pressure at the center uf the re- 
generator) and the AP*s across each of the control volumes in the heater, 
regenerator, and cooler known, the pressures needed in the work calcula- 
tions, P and P , can be calculated as follows: 
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6 




P =P - 


12 

AP. 


v"’’ 


Hoat conduction from hot end to cold end of engine and shuttle loss 
(atep 11 in fig. 3) . - Three separate paths were considered in the calc\ila> 
tion of heat-conch'^tion losses from the hot end to the cold end of the 
mgine: 

(1) Through each of the eight regenerators 

(2) Through the cylinder walls 

(3) Through the walls of the displacer from the hot space to the cold 
space 

The effect of temperature on conductivity was considered. A calculation 
of the conduction loss through the gas inside tho displacA suggested that 
the loss along this path was rmall enough to be negiectraT The dimensions 
used in calculating these conduction losse^plld the locations where the 
te|^perature measurements will be made are shown in figure 9. 

The displacer picks heat from the cylinder at the hot end of its 
stroke and loses heat to the cylinder at the cold end of its stroke. This 
shuttle loss is calculated using the following equation from refer- 
ence 2: 


Q, 


shuttle 


K.7DS^ at 
8CL 


(26) 
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where 

K thermal condactivity of gas 
D dis|dacmr diameter 
S stroke 

AT temperature diifermme across disfdaeer l«igth 
C clearance between displacer and cylinder 
L displacer length 

Conduction ar shuttle calculations are made just once alter steady (9- 
erati<m has been achieved. 

Sum up heat transfers between gas and metal for each componmit 
(step 12 in fig. 3) . - The basic heat into the working space per cycle is 
the sum of the net heat transfer from metal to gas in the heater and 
ei^>ansion*space control volumes over the cycle. The basic heat out of 
the working space per cycle is the sum of the net heat transfer from the 
gas to the metal in cooler and compression -space control volumes per 
cycle. Since it is assumed that there are no losses from ttie r^^erator 
matrix^ the net heat t’ansferred between gas and metal hi the r^^erator 
over ... cycle should be zero. This net heat transfer in the regmierator 
over the cycle ajq^ars to be the most convodoit single critericm for judg- 
ing vdien convergence of regenerator metal temperatures has been achieved. 
However, for small leakage rates, using this criterion to judge convergence 
to steady operation can lead to significant errors. That is, even tlmugh the 
regenerator metal temperatures have converged so that the net heat trans- 
fer in the regenerator per cycle is small, small chaises in the mass dis- 
tribution between working and buffer spaces &om one cycle to the next can 
cause a significant chaise in performance over many cycles. 

Leakage flow between working and buffer spaces (step 13 in fig. 3) . - 
A leakage flow between working and buHer spaces can be calculated accord- 
ing to 
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(27) 


^*^V|Pc 

where 

C constant (0. (K)01 was assumed for samjAe nm) 

’^e'^bulf compression- and buffer -space pressures 

udkmre flow is from compression q>ace to tauffm' ^ace if and 

toe reverse if > P^* This equation was used to investigate the swi- 
sitivity of performance to leakage. This procedure for calculating loss 
due to leidtage irill have to be vqpdated whmi more information about leak- 
age flow becomes available. Relatively small leakage flows will increase 
toe number of cycles required for the mass distribution between the work- 
ii% and buffer spaces to reach steady operating values. 

Compression-, expansion-, and buffer-space volumes (stq[) 14 in 
fllg. 8) . - The cra n ks h a ft angle is defined in the schematic shown in fig- 
ure 10. Compression-, eaqumsion-, and buffer-space volumes are cal- 
culated from the crankshaft angle by toe followii^ set of equations (as de- 
rived from fig. 10): 

Ly = - (e - r cos x)^ s Projection of rod lei^th L on vertical asds. 

Y. = r sin X L_ = Position of displacer yoke 
* y 

Y2 = r sin X - Ly = Position of power -piston yoke 
hr,n.ax=\/l'* 

- ^l,n.ln> + '^e.clearance " E*panslon-space volume 

- Vi%, max 'S’* clearance *Compre8sion-space volume 

Vbu« = ■ '^l,mln> * ''bull, clearance ' BuHer-space volume 

( 28 ) 
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whmre 

e 

r 

X 

\od 


eccentricity 
crank radius 
crank anid* 

displacer cross -sectionai area 

piston cross '•sectionai area minus piston rod cross-sectional area 
displacer rod cross-sectional area 


and and are the two volumes needed to calculate the indicated 
power. 

Work, power, and Potency calculations (step 16 in fig. 3). -The 
Indicated work is calculated according to 




P(dV^ + dV^) 


(29) 


over all but the last cycle. Since heat -exchanger ineffectivenesses, dead 
volumes, and leakage from working space to buffer space are all an in- 
tegral part of the caloilations, the work calculated in equation (29) includes 
the effect of these losses. Over the last cycle, after steady operation has 
been achieved, pressure-drop losses are calculated. To include this loss, 
ti indicated work calculation was revised to 


W = ^ [Pg dV^ + dV^j (30) 

The Indicated power is just the indicated work per cycle times the engine 
frequency. 

The conduction and shuttle transfer losses are calculated only once 
after steady operation has been achieved. These losses increase the heat 
input and heat output by the same amount. Since it is assumed that they do 
not interact with the working space, they affect the efficiency but not the 
power calculations. 
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Th« net heat into the engine per cycle is defined to be the basic heat 
input (hrom step 12) plus conduction and shuttle losses minus (me4\alf 
^e total windage (pressure drop) loss, (tt would proA>ably be more accu« 
rate to use the heater windage loss plus one<4ialf the regenerat- >r wind- 
age loss for windage credit in calculating the net heat ii^imt (as in ref. 2). 
However » if the regenerator loss is considerably larger than the heater 
and cooler losses, It is a reasonable assumption to use one -half the total 
loss for windage credit. ) 

Indicated efficiency is defined to be the indicated work per cycle 
divided by Hie net heat into the engine per cycle (as defined in the pre- 
vious paragraph). 

Convergence scheme for regenerator metal temperatures (st^ 19 
in fig. 8). - The scheme used to correct regenerator metal temper- 
atures between cycles was arrived at through trial aiul error. The cor- 
rection is made as follows: 



N number of iterations per cycle 

AT| weighted average change in gas temperature over cycle for 
r^enerator control volume 

Tj average gas temperature 

I wall temperature 



Tlitnltt 




a. 


wh«r« Fj » 0»4 and F| « 10.0 art tha factors that s©mw to work best. 
Th« lliuii sttp In Ihn oorrtoUon It 


/5ATj^4 

*v,»,n«w *^w»5,old ■” 


AT^ + 3 AT,^ ^ 2 ATg AT|^ 


') 


T -.r /4 ATj ♦ e ATj , « AT, . 4 ATj ♦ > AT,'\ 

^w,e,ntw ^w,6,old J j 

, ... /ATJ♦^AT,OAT,.SAT,.AT,^^ 

*w,7,nnw ^w.7,old ^ ^ j 

^ faAT5MATe>6AT^.8AT8^4ATg>^ 

e.nnw ^w, 8, old ' ^ J -y 

*r . -r ^ ^ ^"^7 ^ ^'^'s ^ ® 

*^w.9,ntw ‘w.O.old *y J — j 


(33) 


In an attompt to improve upon the rate o( convergence, the scheme 
uaed by Urtell (rtf. 7) was tried. This ia a CiU'rection id the form 

'^w.l.new * ^w.I.old ' *‘^1 ^ I * 5,9 (34) 

where Q| ia the net heat transferred from the metal mxle h> the 
gaa node over the previous cycle. An attempt was made, by trial ami 
error, to pick the optimum cmnbinatUm of factors (I - 5.9) to speixl 
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19 eomrergence. This procedure worked in the Lewis model hut not as 
w^ as the previously described i^proach. 


RESULTS GENERATED WITH MODEL 


Some sample performance predictions ami plots of engine variables 
over a cycle are shown for the GPU simulation run characterized in 
table L The assumed metal boundary temperatures are shown in fig* 
ure 11. Performance predictions for this run are shown in table n. 

The ratio of **dead volume*' to the change in total working *space 
vcdume over the cycle is a very important factor in Stirlii^ engine per- 
formance. IMad volume is the volume of the heat-eschai^er components 
plus clearances in the ea^ansion and compression spaces. B; decreases 
tiiie ratio of maximum to minimum pressure Uiat can be achieved over 
the cycle. Better deHnitions of GPU engine dead volume by reHned cal- 
culations and gas displacement measurements, since these sample runs 
were made, indicate a larger dead volume than used in the model for the 
sample runs. Using the larger dead voltime would result in lower pre- 
dicted power than shown here. 

Eaqaansion- and compression -space volumes are shown in figure 12 
as a function of time. Zero time corresponds to minimum compression- 
space volume. The cycle is complete at 0. 02 second. One and one -half 
cycles are shown on this and the following plots. These volumes, t(^ether 
with the plot of total volume in figure 13, are keys to the behavior of the 
engine variables shown in later figures. 

Expansion- and compression -space pressures are plotted in figure 14. 
These plots correlate closely with the plot of total volume in figure 13: 
Minimum pressure corresponds closely to maximum total volume. The 
pressure drop between expansion and compression spaces is a small per- 
centage of the pressure level (<3 percent). 

A plot of this pressure dr <9 as a function of time in figure 15 shows 
that the maximum positive AP (where AP = P - P ) occurs a little 

G C 
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before 0. 008 second. The corresponding flow rates Into and out of the heat 
exchangers are shown In figures 16 to 18. (Additional flow rates calculated 
at the centred -volume interfaces within ttie heat Mchangers are not shown 
here). The maximum positive flow rates Just before 0.008 second corre- 
spemd to the region of the volume plots In figure 12, where the compression - 
space volume is increasing at its maximum rate and the expansion-space 
volume is decreasing at its maximum rate (to give the maximum rate of 
volume displacement ftrom the hot to the cold space). Similarly, the max- 
imum negative flow rates and pressure drop occur Just before 0. 016 sec- 
ond, where die maximum rate of volume displacement from the cold to the 
hot space occurs. Pressure drop is greater for positive than for negative 
flow primarily because of the lower average pressures and, ther^ore, 
higher velocities vdien flow is hrom the hot toward the cold space. Pres- 
sure is lower because positive flow occurs when working -space volume is 
at and near its maximum. 

The differences in inlet and outlet flows for each of the heat -exchanger 
components divide each cycle into two regiems - a region where net mass 
is removed from the component, and a region where net mass is stored in 
the component. These regions are labeled in the plots of flow rates for the 
heater in figure 16. Similar flow -rate plots for the r^enerator and cooler 
are shown in figures 17 and 18. 

Pressure -volume diagrams for the expansion and compression spaces 
are shown in figures 19 and 20. Since the expansimi -space work is positive 
and the compression -space work is negative, the net indicated work per 
cycle is the difference of the areas within these two diagrams (219 joules 
(162 ft-lbf) in this case). The area within the pressure -volume diagram 
of figure 21 represents net indicated work uncorrected for pressure -drop 
loss (231 joules, or 170 ft-lbf). 

Expansion- and compression -space temperatures are shown in fig- 
ures 22 and 23. The temperature variations correlate well with the pres- 
sure variations of figure 14. In both the expansion and compression 
spaces the temperature variations are about 23 percent of the peak tern - 
perature in the space. 
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Plots of indicttted power and efficiency ai'e shown In figures 24 and 2S 
as a functkm of engine speed for three average woritlng -space pressures. 
Other run conditions are as ^lecified in table I, except that for the runs 
at the two lower pressure levels the cooler temperature was 5. 6 K (10^ R) 
lower than specified in table I. 

A number of runs wrre made to compare predicted performance for 
hydrogen and helium. Indicated power for both fluids is shewn as a func- 
tion fre<|uency in figure 26. The power was larger for hydrogen except 
for frequencies below about 17 hertx. (For a series ot runs made with no 
leakage, power was greater for hydrogen for all frequmicies above 5 Hs. ) 

Indicated eftlclmicy for both fluids is shown as a function of frequmicy 
in figure 27. The efficiency for hydrogen was greater than that for helium 
above about 13. 5 hertx, at vdiich point the curves cross over. (For the 
series oi runs made with no leakage, there was no crossover; however, 
die curves did get closer together as the frequency was decreased to 
5 hertx. At 5 hertx, the efficiencies were 35. 7 percent for hydrogen and 
34. 8 percent for helium. ) 

The dida from figures 26 and 27 are replotted in figure 28 in a form 
used by several other authors. The hydrogen and helium curves cross 
over so that below about 12 hertx the model predicts higher efficiencies 
for helium than for hydrogen at a given pressure level. When the runs 
were repeated with no lediiage calculation, the curves did not cross over 
but became quite close at 5 hertx, the lowest frequency considered. Sim- 
ilar curves predicted by Uriel! *s model indicated a crossover of the he- 
lium and hydrogen curves in the 30-to40-hertz range (ref. 7) for no leak- 
age and a significantly different engine. In the same reference, Uriell 
also shows curves for a Philips engine (a heat -pipe -operated swasiqilate 
type) which indicate that, in an efficiency -versus -power plot, the helium 
and hydrogen curves would cross over in the 40- to 4 5 -hertz range. (It 
is not known whether or not the Philips' curves include the effect of a 
leakage calculatimi). 
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Ih figure 28, In the low-frequency range, the efficiency dr<:q)8 off 
rapidly due to the increasing significance of static conduction and shuttle 
heat -transfer losses as the power level drops. 


CONCLUDING REMARKS 

As shown above, the Stirling engine performance model discussed 
herein trades engine cyclical performance. Testing of the ground -power - 
unit mgine now under way at NASA Lewis will provide the basic engine 
performance data necessary for a quantitative evaluation of the modeling 
process for mechanically linked engines. The data resulting from these 
tests and comparisons of the model predictions with the data are subjects 
for future reports. 

There is little information in the literature applicable to the periodic - 
flow processes occurring in Stirling engines. An expanded base on both 
periodic -flow and steady -flow heat -transfer and pressure -drop data for 
Stirling engine heat -exchanger compements is needed to improve the level 
of confidence in predicting engine heat -transfer and pressure-drop 
characteristics. 

It is expected titat better definition of the actual engine dead volumes, 
modifications to the presently incorporated heat -transfer and pressure - 
drop correlations, and, perhaps, tightening of some of the simplifying 
assumptions may be required before the model predicts real -engine per- 
formance with a high degree of quantitative accuracy. However, the 
model 's qualitative ability to predict variations in the state of the working 
gas over a cycle and to predict performance trends has already been use- 
ful in helping to understand operation of the engine, to plan the experi- 
mental program, and to study sensitivity to various engine and working - 
gas parameters. 



APPENIHX A 


DERIVATION OF GAS TEMPERATURE DIFFERENTIAL EQUATION 


The basic gas •volume equations are used in the derivation. For con- 
venioice they are listed i^ain here: 

i(MC-T).hA(r_.T) + C-(»,T, -W-T-)-P^ (1) 


dt ” 


^i 



at) 


PV = MRT 


(3) 


E^^ding the first term of equation (1) gives 

MC- C + r T 2! . hA(T_ - T) + C_{w,T, - w.T.) - P (Al) 

''dt’dt * piioo ^ 

Ufferentiating equatiim (3) gives 

MR^+RT — «P — + V~ (A2) 

dt dt dt dt 

Letting R s Cp - in the first and second terms of equation (A2) and 
solving for 


C T ^ 
^ dt 


yields 
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C„T ^ * M(C - n ) ^ + C_T ^ - P ^ - V ^ 
^ dt P ^ dt P dt dt dt 


(AS) 


Substituting the right side of equation (A3) for the second term of equa- 
tion (Al) yields 


=.bACr,-T)*Cp(-,T,-*„V-P^ 

or 

f - ^ V ^ f = ■ T) Cp(w,T, - w„T„) (A4) 

Using equation (2) to substitute for dM/dt in equation (A4) gives 

MCp S = hACT, - T) ♦ Cp(w,T, - w„T„) - CpT(w, ' *o> V ^ 
or 

MCp ^ = hA(T, - T) + CpWjCr, - T) - CpW^fTo - T) + V 32 (4) 

dt dt 


which is the equation used in the model to solve for gas temperature 
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APPENDIX B 


DrrEGRATlQN OF DECOUPLED TEMPERATURE EQUATIONS 


The equation 


dt 


Idue to 
change in 
pressure 


V dP 
MCpdt 


(5a) 


can be integrated in closed form by solving the equation of state for V/M 
and substituting in equation (5a). 


PV = MRT 1 a 51 

M P 


Substituting 


dT _ RT dP 
dt PCpdt 


. dT , R dP . ‘ dP _ r - 1 dP 


T CpP 




y P 


In T 


t^-At 


= Lllln P 

r 


UAi 




T 


t 
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dt 


WjfTj -T) -WqCTq -T) 


due to 
mixing 


M 


(5b) 


By using numerical integrati)n let 


.pt+At . ^ At 

dt 




M' 


t+At 


f„t+At ^t+At^t ^t+At|,t+At t+AUt+At 

.^t -*1 F , ”i -'^o 


M 


t+At 


At 


M 


t+At 


Since 



have 
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. * (w{^*t{^« - w«*^*T^*^*) At 

„t.At 


rfT 

dt 


due to 

heat transfer 


hA 

MC. 


(T 


w 



J!^dt 

MCp 


(Scl 


Assume is constant over the time Increment (or the purpose of inte- 
grating the left side with respect to time. This is a reasonable assump- 
tion since changes much more slowly than T due to the relatively 
large heat caprclty of the metal. It was a'so assumed that h and M 
were constant over the time Increment to allow Integration of the right 
side of the equation. 


n+.st lUAt 

. ‘ . -Incr^ - TV = t| 


t 


MCp !, 


(T - , 

In — 2 - J]1£L At 

(T^ - T)^ ^^p 


-{\\\ MC )At 

• . fr^ - = CT^ - Tl* e 




-(hA Mt\Ut 
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or 


/ -OiA/BICJAt\ 

This equstioB says that, as the time Increment is made lai-ger, the gas 
temperature iqpproaches die ‘wall temperature asymptotically. Thus us- 
ing large time increments cannot cause instabilities because of excessb 
change in gas temperature. 
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APPENDIX C 

HEAT^TRANSrER-COirnClENT EQUATIONS 
HEATER AND COOLER 


TlM t<|ttaHon usad tn Iht htaler and coolar la 


h,Ao.oom.A*i 

<*i\. 




(Cl) 


which was derived by linearising 


^ -0.0231 

f Dw '\® » 


k ' 


Nbsselt 

Reynolds 

PrwMttl 

number 

number 

number 


where 

h heat transfer coefRcient 

D hydraulic diameter 

k thermal conductivity 

w flow rate 

u viscosity 

A., cross-sectional flow area 

cs 

Cp specifle heat 

This eouatlon is reportedly valid (ref. 14) for turbulent flow of gases when 
0. 7 < Prandtl number < 120, 10 000 v Reynolds number 120 000, and 
AT is moderate. The Prandtl number was essentially constant over the 
temperature range of Interest at 0. 68 to 0. 69. With the assumption that 
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« 0.861 


•qiaaUon (C8) reducM to 


k 


m 



(C3) 


Heater and cooler Reynolds numbers were calculated In the model range 
from 0 to 25 000 but are mostly above 10 000. When eijuation (C3) was 
linear leedy the result was 

— . 001871 ( -2S_ U 15 (C4) 

“ \t**c*/ 


which Is equival«it to e<iuatlon (Cl). Plots of equations (CS) to (C4) are 
compared In figure 29. 


REGENERATOR 


In reference 15, heat -transfer data for a 79 \ 79-wlre/'cm (200 \ 200- 
wlre/Ui. ), 0. 051 -cm - (0. 002 -In. -) diameter wire screen are given. T%.^ 
data were extrapolated to 400 layers with the result shown as a solid curve 
in figure 12. 

A lUiear approUmatlon to this cxirve 


1 ^ 

k 


0.06071 



3.7 


(C5) 


is also shown In figure 30. This equation was used to calculate heat- 
transfer coefficients in the regenerator. 
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EXPANSION SPACE 

This analysis assumss psrtset insulation bstwssn the combustion gas 
and ths txpansion-spaos wall. Heat transfer between metal and gas is a 
combination of radiation and convection. For radiation 

9.ar(T« -T«) 


and 





where 

o Soltsmann ccaistant 

F emlssivity times view factor 

wall temperature 

T gas temperature 

Q rate of heat flow 

A heat ’transfer area 

Vad radiation heat -transfer coefficient 

The overall F is assumed to be 0. 7, 

The convection heat -transfer coefficient is 


(C6) 




Re > 10 000 


(C7a) 



2100 < R« s 10 000 
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"ooiiv ■ 0.0»J(Rrt®-*(Pr)« < i 

wh«rt L U th« maximum distanea from the cylinder head to the displacer, 
and 

^'conv • (Ofeeti number)^* ~ Graeta number > 10; Re s 2100 

(C7b) 

or 

hconv • Btu/hr4t*-®R Graeta number ss 10; Re s 2100 (C7c) 

where Graeta number « Re x Pr x The value in equation ^7c) Is an 

assumed cutoff point (close to the natural convection coefficient). For the 
combined heat -transfer coefficient the values obtained from equations (C6) 
and (C7) are added. 

The volume of the Insulated part of the heater bdoes adjacent to the 
expansion space is lumped with the expansion space. (R is treated as an 
expansion -space clearance volume. ) However, the heat -transfer rates 
between the Insulated tubes and the (las are calculated separately before 
they are combined with the rate for the expansion space. The same equa- 
tions, (C6) and (C7), are used, but in this case the convection strontrly 
dominates the combined value. 

COMPRESSION SPACE 

Since the radiation effect is small in the compression space, only 
convection heat transfer is considered. Equation (C7) is used for the 
calculation. It is assumed that the wall temperature is known. Without 
detailed analysis or test data to identify this wall temperature, it seems 
reasonable to assume that it is about equal to the avera^^e compression - 
.space gas temperature over the cycle. The net result is that very little 
heat transfer takes place in the compression space and the compression - 
space process is essentially adialiatic. 




APPENDIX D 


MOMENTUM EQUATE 

A general form of the conservation>of-mommitum equation for one* 
dimensional flow is 


at 3x 2D|^ ax 


Rate of 
accumulation 
of mmnentum 
per unit volume 


Rate of 
mommituro 
gain by 
convection 
per unit 
volume 


Rate of 

momentum 

gain by 

viscous 

transport 

per unit 

volume 


Rate of 
momentum 
gain by 

pressure force 
per unit 
volume 


(Dl) 


Expanding the first and second terms of equation (Dl) yields 


p ^ ^ V ^ -f V + pv ~ pv^ + ~ a 0 (M) 

at at ax ax ax 

By the continuity equation 


§£+±(pv) = 0 

at ax 


and second and third terms of equation (D2) can be eliminated to yield 
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The first term in equation (DS) is selected in the model. N^lecting 
this term and multiidy the resulting equaticm by dx/p yields 

V av + — 0x + — = 0 (D4) 

®h p 

Note tiiat at zero flow and second and third terms of equation (D3) are 
zero, so that it reduces to 


at ax 

in udiich case the time derivative term is responsible for any pressure 
drop. The significance of this time derivative term could be invest! • 
gated by the use of a comprehensive model such as that of Uriel! 

(ref, 7), 
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APPENDIX E 


PROGRAM LISTING AND DEFINITION OF VARIABLES 

This ^pi>endix incltides 

(1) A listing of the computer program with the short form of the 
printed output for a sample run. 0*he Caleomp plotting subroutines 
are not included in the listii^. ) 

(2) Definitions of the variaUes that are read into the main pro- 
gram using a NAMELIST format 

(3) Definitions of the variables shown in the printed output 
Other key variables used in the program are defined with comment 
statements that ure used throughout the program. 

The program consists of a main program and four subroutines. 

R is programmed in FORTRAN V and has been run (m a UNIVAC 1110. 
All input data are either contained within the program or are read into 
the main program via the 20 NAMELIST parameters. Values of the 
NAMELIST parameters used in the sample run were printed out and 
are snown as part of the output. 
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LISTING OF COMPUTER PROGRAM 


Page is 
CW POOR QUALITY 


Main Pr<^cram 


c P 9069 AH 

C )OCS mill AtlPAItOftl. UPDAieS PISIOM POSXI tOMS— CXPA^&lOH, CO^PffCSSlOM 
C AND DUfrCD SPACC VOiUNCSa INTCGRATtS PAOUTAIVI 10 0C1CR«*INC NONA* 

C OCieNNiNtS NHCN CVCLtS IS CONNUfC* 

OIHCNSION XI2l»VIUI«NeVNOII3l« QNISI*fNAIU 121 

OIHCNSION 101 121 •T6AUll«TNt|3l,fWANGUSI .01111131 
C00N9N OflNCfMUl.NOIlSI.OCtPnSI.Nt^CtNC 

CONNON /CVC/ OeNP.QHCATN.QNCOCN.QCOOtN.OCONP.CNrNTM.eNFMlN. 
lENrCINtrNrNTC.MNllPtNNRCXPt 

2NN«CH»tNOAtCA»NONLCN»UNlll.P»lMIRH.SUNNUNmi.SUMO€Ntl3I.CNlN, 

3NDISP»NPlS1.MtPISItQCXPN«QCXPNtOHeATNtONCAfP« 

OOCOOtN*OCOOtPtOCONPN.QCONPP»ONEONtONCOP.NOlSPPtNOISPN«NPXSTP« 
SNPISlN.OCNONI.QCNONO.OCNOCL.OCNODfOCNOCN. OSHIL •? OCHPA ,TOCHPA 
COHNON #OCUAP/AO«AP. HAI 131 

CONNON /TAVCVC/IOCVC«12lf 16ACVCI13I«1NC1C1I3I 
C0NH0N/C11t/NCUN«»1S£U •1SC12.1SC13 
CONNON /JNDCX/JlP.NOCVC 

NANCLtST /SmNG/P»N£ At6S«riPCVC*t PRINT »I TAPS •C0CFF«TJ,P3.0HtGA, 
PFACTI«rACT2«N0CtC»NSTRI.N0CN0*NU6AS«l6«T0A,IN,RHCrAC»JIP 
DATA AO » AP , AR • APISIR / 

1 6*0t0*0.0*ll .9*00/ 

DATA RCRANK,e.ROOL«N3Cl.01ANO«01AHP,DlARDR.01ANPR/ 
n*S0S*0*B2.1*91»2A*a«it2*TS»2*TS»9*37St0*eTS/ 

9ATA V/3.0 t 2A1 *5008. .033 »S«3 *620 .SAO.PO} .0*3/ 

10 CONTiNUt 

C RCAO lO INITIAL PRLSSURCS. TCNPCRATORCS AND OTMCR PARANtltRS NNICM 
C DCPINC THF NATURC OF THE RUN TO BC NAOC 
RtAO l5.STRLNG«CND=n33l 
0TlNe = U/l0HeCAOFXPCVCI 
NRITC I6.5TRLN81 
C IN1TXALI2ATX0N 

FMULTri. 

IttRsO 
T1NC£0*0 
PSISO.O 
PS10C6:0.0 
5AVCT:C*3 
tPLOfrO 
»RSMO:0*0 
N0ITPt:O 
DO 11 1=1.13 
!FII.rQ*l3l GO 10 12 
TGCVCtlKTGIl I 
12 TGACVCI1»=TGAI1I 

1 1 thctcu I=TNI1 I 
CPCVCtl */<OHCGA«OTlHC I 

C O.CP.CF IN UNITS OF 1N-L9F/ LBN-DCG. R 
GANMArUSOA 
8=9197. 

CP=32557. 

CV:23300* 

8=*32T58 

IT :«NCAS,Nt*8l CO TO IS 
I = *D1M3 



m 


**«»S7* 

It CMIlVWt 
IM|Mt«it 

tmiHtlMf t*!! 

t !*««« «» 
tMtCtMINI 
tUMUftlM 
l«ll«>»»OlVIU 

« 1 1 «K»0«M» *«»(M I «0. }» 1 

DO n |:l«ll 

0Si0S*«St1.4/ |t«t1S*0*OS«OCHt%M 
*i4M«:<IS**/«0«0t 

HtOS/IIOS 

►t*0 

•C*0 

DO to tx|*tS 

to OOIIIxO*«U W «0««t0*«U*0«fORt ttOtll 
ltlOOV:| 
j:t 

IDO COOtUMit 

ir «j , 11 . uiooii 00 10 t»o 
irt It « 4 to,oi,«i 00 VO too 

ir vj.iv.vottH-iv* ootu VO, till 

ISO V 0004 V »»> viw oiioti «<u sv»i vt»,v«»»v'»»,v«,»ii*,%»,*it 
I', tt,» ot DC oom ooouv •,s«,*nti«,t»,»rv»i •• 
tr Vvt.to.uiocvc-in oom vt.il^v 

is< tootiov «* not OIIOI.C on ott ocs oto oit* 

I* ti» on oco oto otio otii ottt 
»»nis •» 

ir t.i,fo,<«ocvct uoivc It, V SSI 

iss roto, V «• not miv ri rj vs vt vt* 

i» Vk *1 V* v« no VII vij«i 
c ti'oitouio or 1000 
1*0 tuo:n(o«i 

in j.ct.iNooc-tii looioviiooim 
ii inioi*i««tti.ov.t,oi 00 10 too 
mir«i 

VVOVixO.O 
DO 1*7 iJtxl.lS 

1*1 «voii.:«vou*«nooi 

tr (jto.ov.oi 00 VO ivo 
ir no .iv, lootovt to ** »f; 

. 11 , itvMvi to VO tvo 

IV I J.tV.ltOCVC'lll DOlVt IttItOI llO{,OUOt«,k,VOI*V.VOI«l, 
ivotnv,vo*iisi*ot,oc,ooiii,ooo«iv,v (tt,Moi 
I to roootv n«,rt.t,rt.i,trt.t«iii,tro.t,ri,a,rt.a«iK«*vo.7,>ro,ti 
IV IJ.tO.VOOOflll ootvc lt«lt}l VlOf ,otiD(o,oivoo 
ito roootv iik,rv.t,ro,i,ttvo,oi 

ir u.to.oocvci ooivt it,itsi iiot,otioto,r 
Its roootv n«,r v.«,ro.i,t<ro*ti 

IV IIVOOt.{0,ll ooiu it.ttoi io,vot,vo 
t»o roootv I* vox*,*i,i}ro.a/,* vtt«*,isvo,o/,* vo:*,isro.o/i 
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ORIGINAL PAOB » 
OP KX)R QUAUTY 


IPP IP50 

IBO ?|ilC:OIIt«EortORVllT(lll 

IMl1|H.CQ*irilC«SPDt|MCI •AHO. J«tP*6l SO TO YOO 
C C*U HCRI CXCNANBCft SUBHOUYtHC fO yPOAYt %k% LUMP MtSSCS* PtOM «AY|S« 
C TCBPCBRTUBCS* BRCSSU4IC ♦ BCBCMCtMOB HCm ICBPCflX YUHrS. 

CALL NCAYX tl 1NC»1C*Y«A •IN •BCVtlO*OII»X tHOOW » J»01 ANO«X 
|B0D«MX,V^SteStB«llM$A%«BNCPA€»BEALBS»C<M4OID»NOCYCI 

t|OttPCsilOtm«l 

C aPBAYt CxBAIISIOM ABO COBMCUtOM SPACE PBESSUtES* UBP/IBPI 
PfOLOm 
pmospc 
POLOXP 

tr U«LT*MIOCVC^SII 10 to 0 %) 

PO|B:IOCLP«A»«OELP96I«0*I*OCLPIY»IOP99ULY^P 

P09Mf:P*l0.So0ELPtTl«0EtPiaimtP(0llArBULY 

PCstOCLPlPI^DCLPnMOCLPtAllorilULY^PPIB 

PC:POOUV<»«OCLPnaMDELPUtl«OCLPtUIIAANUlY 

10 VO t%% 

PAD PltP 
Pcrp 
POIBSP 

Otouvsp 

PAS COBIIIIUC 

ACr|or«PCOLD»/P» 

AC-IPC«PCOLDl/P. 

A:IP^P010I/P» 

DCLPOOSPAIB-PIIOUV 

OCLPNTSPE^PPIN 

delpcl:poout-pc 

C CALCULAVE mVfPBAL IBEOOV OP OOPAlllO SPACE OAS* IPt^LOPI 
UnVAL^O.O 
DO PSO t:i«lS 

PSO aTOVAL::ilOUlACAAVOAlll/lp.«UVOf AL 
PSO COMINUE 

C OUPPCO SPACE BLEED PL0M«PL0T0B ILON/SECI 
IP UPC^PSI .BE* 0*1 BO 10 PVO 
StB)P:*mPrs«PS*PClss*SPIi 
SO VO POO 

PVO SIB)P:COErPS«PC»PSIsp.S/tf 
PSO COBVINUE 

PLOVOOniBSPsy 

IP IPLOVOB.BE.OI EMHMB:CPPVBAII)t*PL010BP0V|HE/lP. 

IP tPLOVOB*LV»Ot EBVHMB:CPsVSSPL0T0BSDVI9IE/IP* 

ENVNSENVN^EMNBB 

SISS:StBS«SlBSPSDVlME 

ABUSI-BBUSI^PLOTOBSOVIIIE 

IP «PLOVOB*LV .O^OI VBAI1VI:TBA||SI*PL0V0B*DV1BCS 
II tS-IBAItSn/MBU5l 
SPO PSOLO^PS 

PSt|)/|BS/IS!(:3ttf«C}'BAwiiitGSI 

AS:|PSOLO»PSl/P* 

C UPOAIE COANA ABBLC* PS t IBADI ANSI » PSiOESf DEBRCESI 
PSIOLDSPSI 

PSI:PSI«0yEBA«0VlNESS»P8SlB 
IP IPSUBE*S.PBS1BI PSt:PSl*S*2SSlB 
PStOEBsPStsSSO.PS.PBSlB 
C UPOAVE PISVOA POSIVIONS IIBI 

AOOVsSQOf IBOOLPSP^IE-BCOANRPCOSIPSI IISSPI 

llOlOrXf II 

RPO.OSXIPI 

Mill sOCBABAPS IBIPSI I «RODV 
VlPi:BCBABNSSUllPSll-POOV 

C UPOAU CXPAIISIOB, CONPBCSSIONt ABO BuPPEB SPACE VOLUNES. Iimi 

violdsviii 

VISOLOSBIISI 



so 


c ctLCiatie «o«K sp*ct cowmi^toN 

C «0««,«*MCNkt TOt»L HOftM, M««H| DlSPlkCC* MODHtHOn^l POMCtt »I$TON 

C MCUURt DROP IIORft lots 10 tROSNSlOH SND tOMROttSION 

C SRSCtSt SkO HRRLCN 

JOLDI/I}. 

Otttt IPlOOkC 10 *OME no 
t* tOMtIO.Ot.O.DI MC«OP:mCIOO*OMEIO 
ir lOHtKO.U.O.OI MUORSHOOil^OltCIO 

HRKCNr .|tRNCNO*OIICNO 
ir I 0 HC 1 O. 0 V. 0 .OI HCMOO:tfCMOO*DHCNO 
tr lOOCMO.lT.a.OI HCNONtltCNOMBIiCNO 

«0«CC<»«H*KICN*I«>«C t«««<III-«ISOCDI»l>* 

OR4MSttRIIM*Ut«(Olt»-»ieiB»*«C«<«<tSI*«tl0L0ll/». 

BM0tSP:«kr«»0>*C**0RI*mi 

tfOISP:N01SP«OttOtt* 

IF <0lt0ISP.«t>0>0l M01SPPrMlSPP*0v0tSP 
IF IOHOISP.U.S.31 OOtSPNSMOtSPMOODISP 
DMPlST:ac•*Ol>•l«l^l•«^OLOtm. 

PPISTFMPISIOIlPlSt 

tr «0MP|St.6T.0.0l OPISTp:MPtStP*D«IPtSI 
IF lOVPISF. 11.0.31 MPISTI«:ilPtStN«OllPIST 
tfFPIST:MFPtSI*«CP«0»P«««)l'K>OI.OI/U* 

DO S» tsi.tl 

irii.ro. Ill GO to tn 

TGCVCIli:tGCtcm*Ti<|l 
til tt«CVCIti:tG«CVC«ll*tG«lll 
tSO INC«Ctli:tNCVC(tl*IN(ll 

IF IPtlOLO.et.J.tPISO .GNO. PSI.lt. S.tPtSO t GO to GtO 

SO to l«0 
c CNO or loop 
SSO COOtlNUt 

C CSlCUlStC SOCOSGC 6«t tENPCMtURCt 0«C« tNC CVCIC FOR CtCM CORtROl VOlUNC 
00 til lit. IS 
IFd.ro.lll GO to SIR 
IOCVCIII:ttCVCItl/CPC«C 
til tS«CVCII»:tt*CVC«ll/CPm 
ill tPCtctiistRc*cit»^For*r 
tCIKPCttSCTCUl 
tGCPP«:tG*CVCIIll 
00 str tiFi.G 
jj:||*s 
tN«ill||:0.0 
GGt tNSM Jji:tOOO. 

IF <J|P.6T.O.«NO.J.NC.<NOCTC-III 60 tO 660 

IF < J.lt.lNOCVC'in MRIU <6.1611 tlMC.PSIOCG.K.P.Pl.tGkU I. 
lt6Mlll.Pt,PC.PR|M,PR0UT,« tll.VUll 
IF IJtP.CO.OI MRIU I6.661I IS.tGt.tN 
660 CORtIMUC 

IF< J.lt.MttRt .OR. J.Gt.NOCNOI GO tO 6tS 
C CORRCRSCRCC tCNCNC FOR RCGCRCRRtOR NCtRl tCRPCRRtURCS. 

00 6t0 JiFt.V 

6F0 OtR<Jl|sFGCt|«OtNIJtl«FRCtmUHNUI«ltlFtUNOCN<Jll 

tHttl:tN<tt-<S.«OtN«tl*6.60tN<6l«S.*OtNttl*2.6DtN«OI*OtN|RII#S. 
IN<6l:tNI6l-l6.«0ttl<tl*0.6OtN<6l*6.P0tNttl*6.PDtM0l*l.*0tlM6ll/l. 
tNlti:tP<tl-t|.6DtN<tl*Z.60tN<6t*S.60tN<ri*l.POtN<GI*OtNI«tl 
tN<6l:tN<0l-IF.P0tN<tl*6.6OtN<6l*6.60tNiri*G.POtNIGI*6.P0tN<0l 1/1. 
fM<.9|:tN<6t‘(01P<tl*l.P0IN<6l*S.60tH( tl*6.60tH«GI*S.«0tN(9U/l. 

GO to 6tt 
6tt CORtlMUt 
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ORIGINAL PAGE IE 
«)F POOR QUAUmH 


00 OT) l:l«S 
07) aiMltsO.O 

tV5U&P:o*SUM/rtO*1 1 no 1 7 PC I 
C*Ul C7CI niHt«UT0m.$7P0«i».*»6MSP» 
00 0)7 ):itU 
IF Il.tO.tSI 00 70 0)0 
70CVCm:0.0 
S)0 76kC7CIIK0.0 
07? 7NC7CU)O.0 
PROUN70.0 
• 0|7PC:3 
j:j*t 

If 7 J.«lC.7«IOC7C-in 00 70 080 

IPL0770 

I7C8S0 

0*VC7:7|Nt 

iiocsn. 

080 r «J.LC.>IOCVCI so 70 120 
730 C9H7INUC 
SO 7.. 10 
1300 eC7T»XUt 
OTOP 
E«I0 


Subroutine HEATEX 


C 7MI0 0U880U71NC CALCilLSlCO HC*7 lOONorCR AND UPQA7C0 7CN8rOAIU»CO ANO 
C •AtOOUAf 

OU8KOU7IHt HCA7X 17 INC .78»70A.7NA»RCTM0.0RtA»A007« JC«CLC,01AM0, 
mNlN.ROOFNXtVtOieStV t*7HSAO<l7NCPAC»RC AL60»C0N0I0tN0C7CI 
COANON 071NC»r(l2l,8BU5l.0Ct.Rtlll,P«RC»RC 

CONNOR /CVCF OCIRtONCA7R,ORCSCN»OCOOLRtOCONR,CNrRTN,rNPNTR« 
ICRPC7A.CNrATC>NRRRtttRRC«P» 

28RACNP,URAlCX»URAlCN,URRLP.HRRH,0UNNUNi:SI ,0UN0CN7 1 5 1 ,CN7 N , 
)ttD10PtUP107.N7P107tQeiPN»OCKPP»ONtA7N,OHCA7. , 

A3COOLN,QCOO(.P,QCONRR,OCONPP,QRC6N,aRCOP»UOlOPP .«t010»N »ttP107R. 
0«P|07R,OCNORI *000080, OCNOCl.OCROO.QCNOCN* 00N7l.*7SeiPA*76CNPA 
CONNOR /0C17AP/A0*AP, NAIDl 

OlNCROIOR COCPAINOl.OPAINOI 

OINCROIOR 76AI1)I, 7NAtl}l*76A0ll)l*7NA0n)l» 1121 

OlNCROIOR CPMIISI, emi,7Nl«l)St,767l2l*OM«ni. ACOI))). 

lANIIDI ,COROIl)l,V10llSI*RC7lllO7l)7» OR IS I 

OlNCROIOR Vllll,MOOLOI)'ll*HII)l.rAV6l))l*OR077:i)l*rRIC7l|)l 
3A7A DN/I.NST,)*0.II90*SN0.00S800,SP3.0N0«,1.N87/, 

I ACO/).?00,)NO,0006I.S*O.N)80,)NO.OSI2N,).ISO/, 

I AHI/).3*Zno.72'« .. .021, 0*82. 0S*S«2. 288, 0*0/* 

I CPN/N*)0.0.0«.0ll. ROIO.O/, 

I U/0, 0,20). 08, 2. 008, 0*0. I?8,)*0.! 87,0.0/, 

1 IOEI/1/ 

C )NI1I:N70RAU11C OlANC7{R,lN.i ACSIlIzFLOM ARCS, 1R2| AH7INCAI IRANSFCR.tNll 
C CPMIICMtlAl NCA1 CAPAC177, 87 U/L8N-0C6.R t XL I II :lCRS 7NS roR CAIC '.'I.A71N6 
C PRCOOUAl DROP, IN 

C CP.CX ARC IR B7U/I.BN-DC6.RtSA$ CORS7AN7, R, IS 
C IN |R'tBP/U8N>0C6 R» 

B:. 02)08 
CPtS.NBN 
6 anna:i.)9n 
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9 z 919?0 

ir 60 TO S 

SS.01613 

CPri.^39 

Or«63?« 

6A*INAr|.6^7 

5 CONTINUE 
£v:CP/SANMA 
00 13 i:i«13 

IF INUGAS.CO.AI 60 TO 6 

C CONOUCTltfllVv CONOITI NAS UNITS 8TU/IN»StC*0C6«0 
C VlSCOSITVtVlSm NAS UNITS LBN/lN-SeC 
C090II»S2.Aeit-9»T6Alli«|*?63£*6 
tflSlI l=«.6e3£-l09T6AUM2*6l9t*7 
60 TO 0 

6 C ono I n : 1 •96^c -9#T 6A n I • 1 1 SE -6 
VlSlti:l.0SeC*9*T6At]MS«919:*7 

8 TNAOlIt«TNAIII 
10 TBAOIIISTGAIII 
P0L02P 

C CAICUIATC P0ESSU0E«L8F/1N7 
SU9UT :3«3 


SU«fV$0«0 
00 1« t:l»13 

SU9Wl=SUHUT«U6iIIPT6All I 
19 SU9V:SUNV«VI1 I 

P:8tSUNMr/lSUNV-B#ft9U«ll»-SI63l9frEAt6SI 
F6AH«|A = tP/P0L0IP«l ISANNA*! .I/GANNAI 


00 ?3 Jrl,l3 
T6A0I JlrTGAOIJlPFGANNA 
?U TNIM JirTGAOl Jl 

C calculate HASSES, UGCII, IN EACN GAS LUHP, «L0H 


SUH:3«0 
00 700 l:i,13 

TOO SUH:SUN«PPV| n/ldPlTNlXflMBAP TCAL6SII 
00 ZZO Ul,13 
220 UGOtOtTiruGtl » 

00 790 Irl,l3 

790 U6UI:B9C 1* -S16319P9VII |/IP*SUN9l INIX I It^BPPH^EALGSl I 
IF IHOEX.C0.2I 60 TO 793 
00 797 1=1,13 

792 UGOLOUIrtlGU 1 

CALCULATE FLOU RATES, Ff 1 1 ^PLBM/SEC, BETmECN GAS LUHPl 

793 HOCX=? Vi 

F « 1 1 : f UGOtO 1 1 1 *U61 1 1 1 /O TINE 

30 300 1=?,17 

300 FIIISIUGOLOII I*UGI1II/0T1NC*FI|HI 
00 29 1=1,12 
tF|Fltl«LT,0*l 60 TO 27 
TGClirTGAOIll 

IFC1«6T»9 .and* 1.LT«10I TGI 1 l = T6A0l 11 -I THAI S I -TNA |9 » |/ y 
IFI t |#DT!NE9ITNAISl-TNAl9n/6»/U60LOll» 


GO TO 79 

^7 TGI ll=TGAOII«li 

IFII.GT.I •AND. 1,LT.9I T6 11 iriGAOl f *1 1*1 THA I S I «^THA f 9 1 »/8, ♦ 
;FIII«DTIHC9I1NAISI*TNAI9>I /e./M60L0ll*n 
79 CONTINUE 

C CALCULATE AVERAGE FLOU RATES FOT EACH 6AS LUNP. 

FAV6llt=Flll 
00 30 1=7,17 

30 FAV6III=IFI1-1MFI1 l»/7. 

FAV6I13I=FI13I 

C REVISE GAS TCHPCRATURES TO ACCOUNT FOR NIXING* 

IF |FI1I*LT.0*0I TNlXflKluG0L0lll9TN|Xlll*Flll90TTNC9 
ITGini/UGIll 
00 100 1=7,17 
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_r ., 

v)F POOR QUALTW 


100 fNixiiiriiisoio«ii«iMtiiiii«irit»ii*fsii*it-rni*t6iiii« 
;oMNci/ootii 

ir iru?|»ST«0«0l ?Wtlltt)iriMSOt0tl5l#TI»lHIISI«FCl?l# 

i9iiNc*fGim»#oeui) 

110 coxfimif 

C COLCULOU orvootos no S» oeviiottls NCOf lOAIlSrCR ARC AS* AMY III; 

C 4CAY YQA^srCA COCfF 1C1C«TS» NAIll BCYMCCM SAS • HCYAt lUMPS. 

C MAIII HAS OMIYS BTumC*0C6«0 

BCVH0ni:0NlllPABS«Mlll/IVlSUlPACSUll«»333l 
AHf Cll:S«lAl6P0IAH0PCXIll-VlNtlll#l6.b 

HAI’ |rcotlO«ll/ONIllPCt«BYir-SPOMI||PABSIMll I/IBISIII# 

* IIMIS.IPANYCII 

C YO N : CXPAHStOH SPACC PROCESS AOlABAftC SCY MB 1 :0*0 
C HAIli: 3.3 

C 

C QUEHCHIIIB IH CRP* SPACE— RADI A Y1 OH .70 EFFECYIVE— H€ONV = S/3bOO/ lAA Ml 
tri 

Rt(li:Mll*AlNlll«.01 

KLI1SI=VI1SI/AP«*31 

F0A£0*7 

irCABSlTMIXIll*VNAIltl.iE.1.31 BO TO 02 
3RA0:0.17TPF0APMT1IIRU 1/100* YPPO-I THAI II /lOO.IPPOl 
HRAO=ORAO/ 1 TM IR 1 1 1 -TMA till /S63 3 */ 10% * 

BO TO 01 
02 NRAD^0*0 
0) COHTIHUE 

PRZCPOVI Sll I/COIIOU I 

HC0H«r0*02SO RETHOl t lOPO* BOPROO0*O#C0ll0l II /ONI 1 1 

triREVIiOUI*61*?100*0*ANO*REVIIOItl*LT«10000*l HCOOtf : nCOM«p II * « CONI 
1II/RLIIIIOpO* 7I 
BRACT2=REYHO(tlPPRPONIII/RLtII 

lPlREytlO«ll*LE*2133*3*A«IO*BRArT2*BT*10*l HCOMRzl • 06PGHAET2 #p3 .ITSp 
ICOODIII/DHIII 

lFlREVH0<tl*LE*2133*3*AN0*GRAEY2*LE*10*l NCOt|RrS.3/3633*/10« • 
tril*E0*13l GO TO 00 

C ACCOUNt FOR NEAT TRAHSFER TO INSUIATEO PORTION OF HEATER TUBES 
C IAOJACEnT To erpansion SPACEI along HITH neat transfer to CRPANSION 
C SPACE NALL 

«Ll:l.3?l 

REVN0i:0Nl2lPABSIFIlll/e«0/IVTSIllPACSI2IM*000l 

ANT|:10.9SSPRL1 

HCONVUO*02SPREVNOlPPO«BPPRPPO«OPCONOt 1I/0NC2I 

HAUINRAO«HCOMVI lOANTI 

HAI 1I:INRAD«NC0NVIPAhT1|I«NA| 

c 

00 03 1=2*12 

RETNOf II^ONf I IPABSIFI I-lMFllll/2*/a*/IVlSI 1 IPACSIII I »*000l 
IF ll•LT•S•0R•l*BT•9l NAI 1 iSCOftOI I l/ONII l#l | .e7tC -S p 
lONi ^IPABSUFI 1*1 MFf 111/2* l/B*/ 

7 i VI SI T I PACS 1 1 1 1 « IS* IPANT 1 1 IPB* 

IF ll*GE*S*AN0*t*LE*9l HAT 1 l=CONOl II/ONt 1 |Pl3*0637tP 
lONI 1 1 PARS I IF 1 1-1 MF IT 11/2* l/B*/ 

71 VISI ItPACSI 1 1 I♦3*7IPAHTII |PB* 

00 CONTINUE 
NAIOI=0*3 
00 N1 1=S*9 
N1 HAlIirRHCFACPNAIll 

RETNOI l3t;0NI 13tPABS|yil2l I/IVISI UIpACSI I SI I • .033 1 
ANTI13I=G.?837pDIANDPIRODVNR*ROOVI«S*S 
HAI13irC0N0ll3l/0Nll3IPIl*a7ir-3P0HC 13IPABSI FI 121 1/ 

1 1 VI SI 13 IPACS 1 131 MIS* IPANT II 31 
C TO mare COMPRESSION SPACE PROCESS ADIABATIC SET HA 13 =3.0 
C HAIlSir 0*3 
1=13 

tri|*E0*13l GO TO 02 
NO HAI ISirHCONVPANTIISI 
C 



H 


c cticiitftit •Air% Of Mcof loinsftft ocitfccti SOS ««cfoi our n«m 

C )AS to OtCOaM foil HCOI tooiisrcil to 00 ftoo Mfioi* 

00 110 021,13 

MOtfCS-OlINfONOlOI/lllOIOlOCOl 

fCItsU^tiOINOOtt 

010*21 toAOlo I •toil 101 lore fo 

01012 - 010 * 000 loioco/orioc 

orcttrrciooooioi 

0010*2010*000101 

lUOOOOl 01 21000001111 «00 10* 

*UOOtOtOl 2 $oOO€OIIII#OfC?* 

1?0 IO*IOi:lO|MOI«OtO* 

If iitmf.it .monrc-sii oo to oo 

C COtCOiOtt Ofosttv.rotcf 100 fACtOtS, AOO ofit* OOCSOOOfS 

c roo coco iuoo« 

00 ISO 121,11 

ISO OOOItltirO/lOOltOOttMOOfoOCAiOSt » 

00 100 |2?,u 

If ioevooiti.ii.ioo).*oo.i.ii.o.9o»i.ot.oi fticfiiitio.m vooiit 
If lofsooiii.of •isoo..*oo.t.it.s.oo.i.of.fi r»iciiti 2 <i.ooo/ 
IQCVOOl liooo.i 

If ll.OC.O.OOO.I.iC.OI r 0 ICtltl 2 S. 00 ock 0 t- 0 . 0100 oOCni 01 ttM 0 .S 0 

100 OCiM ti 2 roi€t iitor*«tit io*ostr*ottiiiokiiti/ 
it SS.ioOMt t looostf 1 1 lOACSt t tooM 0 »O 00 . 1 

c mcssuoc 0000 c*iCui*tiotis 
00 101 |21»00 

101 OOktt t^OAOO* 

*CUkO:ot lir«0O|.l0l0O7.?S/0. 

*CSCOO:ot lSI/*ros.lO|OOS.7S/0. 

i2| 

irtrtii.if.o.oi its 

fro 

roirrtn 

»Sf2Mft 
roo^f toi 
f902M0l 

ril:rttfi 

CALi «OCLMi»*CSCkP, ACStO f MOOtSfl »U tOf » •OOSt t Of t ,0t StSf I troi 

ktf«corrirtoiitOPktoii,ii 

c 

CAti «0fiOtS#*CStDf I•*cst0fl,0ltt^^l•kit0f I » 00 StVt 0 f t,tistof itroi 

««f»corMtofi*oPkiOft*it 

CAU kOriPtO.AtSlOf l•A€StOfl•DNt^f t*li tOf l•0•lSttt0Sl•VlSt0St,r0f 

«tf»cocritosi»op«tosif St 
c 

OPfOierO.O 
iif 2hiifi 
00 lOf |2f,lf 

tn i.ro.n «itf irns «i.sf i 

irtt.it.S.OII.I.St.OI OttPCriO 
ini.Or.S.AIiO.l.iC.OI MVPCril 

CAll kOCtOIAtVPC,*CSI|lt*CSttltOHf ll,ntl 1,003 tv<tl,vis 111 «r Avon 
X,f,COfflt l^f l•OPItt«^l»lMfll 
kitftrns 
Sfirtli 20 o«n*ii 
oi^rtiC20oroic*DCiPni 
los cootiour 
c 

trtrtol.At.o.oi 00 to *i 

CAii koriot I. Acstooi t.offooOtPNisoi ,Ri tool ,D03tvf not ,vi\tooi ,roo 
itftcocrkt isi, jMtisi* isi 

*0 to Of 

01 cootiour 

c*ii tor LP 1 1, .Of. ODD, AC St 101 tON 1 1 01 fkitlOl ,003 1 ft tonvist loi.roo 
x,f*cotrRt I6i,i)rvmi,i*i 
Of cootiour 
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C«U K0rt^lS»A€S«U»»«0%SIVO»>Nll2KHUin«0NSlVM7l»Vl$n?»tr t? 

«« i*torrii« IT 1*01^111 1 ?>•! ft 

CAlt mi»tJ»ACStlll»ftCSCON« DMtIf l•VISlm«r 1? 

i»i»cocr«it 9 t*tmiiit 9 i«ut 

mutl70«0 


00 OS t7l«|f 
OPSM'« 7 a»sy«t«o^s u I 
01 COOflUUC 

OCi^iiKo^rtit 

DCcrtmrOOSUR 




0 rL»«C?l 70 (l^l 0 ?M 0 OltQ|M 0 ^&( 0 ft 

OCl>^tQSl7D(i^tOSMOO«tOSt 

CCIM 00 1 :OCt^ 130 1 1 1 SI 

OCL^t lOirntL^nOMOflttOI 

OCL^I l?l70Ct^tl?M0MllTMO^SIltt 

I TO COOltOUC 

C CiiCUlOfC HCAf COHOUCtlON iCOit HfOt CONOUCTtOO SU600Uft«irt 

tr tOCOOO*lll. II Citt CMOCf ttM0«0COC6 •OCn*OCPtlOO»OC«N*QSNf u« 
MOOASI 
OCOOOtl 
00 COOllOUC 

C OCOlSt OCOdllOOYOIt MtliL Uf«nit*lUO(S. 

00 If? jrs^o 

Iff tnoi jtrfoit ji«giai*Dt mCfc^t Jt 

c coNvcoi oijMi 10 rf^tor i sun uo roo cocn cooroocoi* 

0t0l70«0 

OC 173.0 

SMCSttn.Q 

OHCTOrO.O 

0HCIO7O,3 

OOCOCtO.O 

000070*0 

000070*0 

0C00t:O*O 

0C0lOr0*0 

OCOlO:o*3 

OCOOrO*0 

00 Iff J7U1S 

OODD7QIJ|ODTIMtOf ff*s 

tr ij*co*i> oiK7ack«ooon 

tr tj*or *?*ooo*j»it*oi 0Mr«f70NtoT«Qooo 

ir tJ*O(*?*«OD*J*U*0*0O0*0000*tt*0*0l 0NCfO7QN€IO«Q«00 

tr « J*0( •^.OOf)*J*iC*0*OfiO*0000*01*0*3V OHCfO JHCtO^OOOO 

tr 4 J*Gf *s*000*«l*ic*0l OOCO€700(Or«QOOO 

tr 4 J*GC*S.AOO*J*WC*0*OOO.OODD*01*3*3I OOG^70000*OiDO 

tr 4 J*Or*^*OOD* J*ir*0*000*OAOO»it*0*3l OOG07QOGO«OAOO 

tr 4 l*GC*IO*AOO»J.U*l?l 0€OOi7PCOOi»OODO 

tr 4y*Gr*n*OMO*J*tC*t2*000*0000*LI*3.0l OCOlO^OCOIO^QADO 

tr U*0r*lQ*0O0*J«U*lf*Oli0*aO0O*Gt*0«OI OCOiO:QCOI»»OOD{) 

tr 4j.co*iii ocooiocoN«o*oo 

Ifo 9f0l70?0?«0000 
00 Its J7|,s 

Its 004 JirQOt JMOIJlODttOC/COfMJI 
0(«r:oCKt«QCS 

ir lorx *L?*3*0I OCXOM:QtlOt)i*or.X 
tr 40CR •6I*0*0I QCXOP:QCXOO*0(X 
OHCAtOrOMrO?««OHCO? 

ONIAtOrOHroro^QNCtO 
ONCOfOiOMCOtr^ONCtO 
ooc Gl 0 700r Gt 0 «00t GC 
90CG07QO|G«I«OOGN 
OOCOOrQOt 00*0000 
9COOtOiOCOOtt*OtOOi 
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Q c*oi « Item ti ••mu 

0C0<II»?««0««9 

tr ««C*N>tt*0*0» •CONMSOCtM'MOCMt 
aC<lO*|:OC<lO*l ««€•!« *0ttNC«»VV>S 

•c«o«ox«c*) 0 *o*«mt *otiw*fvr*i 

•c«oci t •c«oci *«c n«ot INC* t n. i 
•cn0o:«c«ioo««coMeo*ot iNC*f n. I 
•C«OC« :«C«OC<« *«r «M«Ot |MC**?Y.l 

no coottiiMC 

mr*tN:ciirotN>««UNiir<«i*co,)*>i«to««»*ottHC 

too 


Subrouttne CYCL 


c tun suooootitic i» c*ivto j«%» o*tt ho ih c«cit n cooHt't. 

tit cotcuiom— OCT HO* 10 000 o*»». on noon ho c»Cir» looicouo nmeos 
c OHO trrtcitocv. nc. otot ooo nooo auoottttct oot ttooro roo memottoo 
c or aorottn ovto % rtai%. aruo >tooioo, tortooatto ocat aoo oooa 
c ouaotitiis ao( ontt ro tcoo to ooioaoartoo roo otai rrcit. 
tuooouvtot era ttiot .utoiai.stooav.aoonaoi 
ototottoo «aoao«i,stoonii.aai 
cooooo oiiot.r 4Ui,«aii$i«Mio«lsi,o,0|,oc 
Data j,«,$Toor.Oi.oiio/r«o,o*t*a./ 
aata jtotror 

cooaoo rcrc/ «no,oi«taio,oo(aro.acoowo,ocooo,tororN,torMto, 
Uorcto.rororc.ooar.MOacao, 

roaacoo,aaaua,«oaiCN,aoaio,M««N,twooMO«ni,svNOfo«ilt*l otN, 

laotif «ooia,«ioui,o(«oo*oc«oa,o<«taio.OHato« 

*OC» 0 (.<l,OCOOiO,OCONON.OCex?<',OOCSO.OOIOO,vOltOO,HOISOIi,« 0 |Ma* 

saomN,ocoooi«oroooe*ocMaci«ocNOo.oCNOco, osoti,totioa«tocooa 
cooooo /ta»c«c/ioc*cm»*t«acfC<i>i%tNcvc<in 
cooooo /joota/j|a,)iocH 
fowi«ait«ct «aaotii,of «oi 

rsa*t 

JI 0 |:j| 0 |*| 

ir«oi.t>no.»t.iiac» oiono:o. 

oaf|o:t|MC*OiOt|M 

•l.Ot|o:tto| 

0 to: 'Otar -OMc a to*ociiooi •acooa •Qcooe*ocooco*aaMU -«ooi t h. 
Olot:*ocao-OM(ava 

aout:acooio«ecooo*ocoooo*orooCi*ocMio*ocooco*oiNii 

ooaL I tMOai t a*«oairn 

««aaaerMOaTOl*tioait 

trriot:MOatot/olo 

(rra:«aaa/aio 

tf rr I :f orotM/torNio 

ttrreuorrto/iMroic 

•«tMa;woaM/f%s3.t|)(ii|i«i 
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PM«Kt»:PN»HP*.T«S7 
meosi./ocLTiH 
eMSMssr»0il»*FftC«*69* 

c stene ttRuaics res catcc* aueN of atcRaars oacR s cycles. 

SIRRCla.llsRIN 
SfORCU.tlSROUf 
ST»RC<R*SISMR«e>F 
srRREiRtaisRRaCHF 
SfRRtlatSiruRMVRf 
STRRCIR.aireFFfOT 
SYRREta.flrRCFFl 
STBREIR.ttsRCFFE 
ST9RCia.»»:9aCKIl 
Sf9ReiHtl3t=FR*NP 
simia.msFRCo 
ST9REIa,l?|:aRKLt 
St9ReiR.lS»:CMlN 

STBRClK.taistlYOTaL 

ST9RC«a»lSi:CllSPH 

SfORC«Rtiai»«Si>9R 
StSRCIR.iriSMOISF 
sfmia.iaisRPisi 
Sf9RCta,l«|:MTRlST 
ST9RCIR,90I:HRRLF 

StBRCla.tllsKaPN 
STOREIRtZZiaOCRRP 
ST9RC«R.2Si:RNCaUl 
Sf9RC«a,2«i:0NCaTP 
St9REta.2Si:8COOLM 
ST9REta»2ai:OCOOLP 
ST9RCta,2TIS0CaMFN 
SI8RE<a»2ai:OCOnPP 
STgREtR.29|:BRE6R 
ST9RCIR,30I=0RC6P 
STBRCIR.SIISIIOISPP 
ST9Rei«.S2l=«01SPN 
STBRCtK.JSISUPISIP 
STBRCia.SaiStfPtSUl 
STORE ta, SSI SBCNORI 
ST9REtH.Sai£QCHORO 
STBREia.STlsOCNOCL 
STORE IK.SaisQCNOD 
STBRCtR.9«ISQCNOCN 
STBREIR.aOISTSElPa 
STBREta.RlIsTeCPPa 
ST9REIR.a2lrOSHTL 
STaRCIR.RSIsOim 
STOREia.MIStfRHBaS 
BO 99 Lti.lS 
tni.EO.lSl 60 TO 91 
R:L*9B 

STOREIRtRisTGCVCUl 
91 r:l*60 

STOREIR.NIsTOacvCILI 

nn*?s 

90 STOREIR.NISTRCTCUI 
190 CORTINUE 

IF I JtP.OT.O.aMO.JIPl.NE.NOCTCl SO TO SOI 
IF IJIP.6T.0I ORITE I6»202l 
202 FORNaT <• LOST CYCLE* I 

YRITE 16,2001 OtN,OOUT,URaERP,VR«CNP.HRaTOT,CrFTOT ,9CFF1,REFF2, 
•BRE6 ER,PmRHP,FRE 0,URNLT, rNTH,UTOTai, 

•ERSPN,aV6MSP, tfOISP.MPlST.OTPtST.HRaLP. 

•BERPN,OEaPP,ONEaTN,OMCaTPtOCOOLN,OCOOLP,OCOMPR,OCONPP. 
•BRE69,0RE6P.H01SPP,M0ISPN,IIPISTP,llPtSTN,0CN0RI.0CN0R0, 
•9C90CL,OCNOO,OCNOCN,T6ElPa,T6CHPa,OSNlL,OlNB,MRaBaS 


ORIGINAL 

POOR (jUAIJTY 
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130 roMtt »/• OI«is*,r».l,l»»*OOU»S»,rO,l,|li,»ll««C«^s*»»*«J»l»» 

j»tmo»s»,ro.i,u,*wrris»,ro.i.ii»»ofrrit*,ro.s,i«^ 

S* OOItCOS*»rO>Si* roco:*trB*s* 

•• «o«iT:*,ro.>. • cotHs».ri.s.» otoiois*, 

c«n»os**rio.».* *»oov»s*»r«.i* 

*• ooij»s**r6.i,» ooi»»t».ro.j*» o?out$*«ro.s,* !i»«iO:**ro.}/, 

T* ot«»os»,ro,j,» •MC«YNs**ro>s»* •MC»f»5»,ro.j, 

•• ocooi>is*t*>*>** ocooios*»fo«i»* oc#ooos**ro«i»* ocow»ps»,r#«j/, 

f » ooM ^ s *. »•.!»• ooisops '. ro . j .* 

I * MU $ fos ocii *# is ** r *» J #* oewoos * ♦»»•»»/ 

»• OC«OCLS*.M.I»* OCllOOi»*r$,l,* OCNOCHS'tro.S.* »Ot«Ms*,f»,J, 
j» otMUs'.ro.i/** otMO s»,ro.s*» momo«s s»»ro.J* 

SOI CMtlNVf 

00 soa isi.oo 

S30 ««*IIi:9.3 

tMJ.rO*S) M to *00 

•CfUON 
«00 «:0 
j :0 

00 *S0 ««:|»00 
StOOtlll»«M:SIOOC«ltRRt#0* 

00 *00 

«tO SIOOini*MKStOOCIll*0«l*STOOt«t>«ORI/S. 

SIS coMtiNur 

tr <jio.or.o.oHO*Jioi.«r>Nocvei oo to ooi 
001 ti lo.oooi 

MOttC U.rool lStOOCIll»R«l»iiN:l*«OI 

ir (jio.tt.ot Mom lo.osoi 

30 * rOONOT I* OOEOOOC TtOOCOOtUOCS e«CR LOST cvuc*t 
oottcio.toii focvctiooctc.incvc 

301 too***? t» t»C»Cs».*«»l3»0.0/.» tO*C»CS»,lSfO.O/.» tl*cns».I wo.o 

too roooot iiNO**ovtoast voiucs or lasi $ cvcics*! 

031 OttURM 
CRO 


Submitiiie CNDCT 


lU9t»0UYt«ir CHOCT ltl8«»0CAC6 tOC Vt *0C08ID0*0CAI|»0SMYU »8MI6AS • 

fwau ji» irfim*fca«M 5 i% 

I acaiii?it OCAMi^i 

CC«IN 0 M /CYC/ 9 t*^*OH|atl»tO*>l 6 CII»OCOOt«*OCO*lPttl 8 r»fM»t«irMfll» 
it^rcf 9 *t«irQic 

s « oisp * 9 pm *» i ^ nifOri ^ K . 9 C »^^* OH € aTii »« Meaip , 

a 9 COOC* 8 *OCOOi^%QCOHPl 8 , 9 CO«P^fOW 68 l* 0 »lS^* 90 IS^^tyOlS^WtliP|Sf^* 
SyPlSYytOCypPl •9CllOIIOtOC8lDCl*9CllOO * 001001 # OSMlL#ISC«Pa»TOCHPa 
OiTa ioD#tai*to 7 #?cYi/i 7 ao*#im##»Ps##i?ao.#ii»o##Pt$,/# 

I ai/,a;,a j*oa,ai j,«ta/i*jo#i* 6 *#i» 06 #u§J# ^•lPo•o•aoo/♦ 

1 aii,ap,ai#ni#iii?#/o.asa#o*o?i#o.oos#o#aoo#o*oto/# 

I acoaoo # ocoi 800 / 0 « 50 / o # i # n 6 /# 

I TcaN/iaoo»#iia 3 * tiioa*/# 

I acao*Dcai«/ioiS 9 »t*lftt>.s#o*^&^t 

I CLcao , onPo » stooae /# oio « p * Y »« i * n / 

c acoNOo*acYt#aco«ioo#acay*-*aoca % yn 1 N^ 
c 0CoiiOR#0CYit0C08i0Of0ca8i-**0iita8iCf i» iw 



ORIGIN M' 
OF KX)R 


c tcvk.*tcM»- iCNPcaiiuRCS tN otsms * 
c CUM— cic«*«*ice ktmcch oisnacc* mo cviinot*. in* 

C OIMO OlSOLACtt OtftHCUO* |M 

C SIOMt OISHACCO S100MC* IN. 

C CNtCULNIC MC*T CONOUCICO YHR0U6M CVllMOCN*BYU/»CC 
TCN«ei=<tcnui*TcvLUiwx. 

ICANS7SI icnui *icvu II *n . 

lHCNOi:i.OO>SY*tC*»tl*«.**$l/IS4aO.*U.I 

ir I1CNVOI.OT.IBM.I TNCN0t:U3BN4TNlCM»t*S.*Stl/OM0.4II.I 
TMCNO?:«.OOieT*TCN«B>*4. 10)1/1 1400. NIB. I 

ir MC4VBB.BI.1040.I tHCN0B:i.)9«4Y«YC««BB*t*«SBI/l )400.*IB. I 
BC«ils4.tBS«lMCN0IMI>*<BSHIBI4l«CVMII-tCVil>ll/l>iS*4i0SIIIR)- 
l•I>l•IRB«*t^ll/lt•S*«IBI•ll)t-•lXIIII 
IC«l.t:4.XBS4lNCIl0Z*BIX*IB«-l»l4IIC«llXI«ICyLI)ll/<nN««l04IMt<i> 
IBItlNIIDMtBII/IIIINNBIBiyiBa.BIBIIII 
BCVt:iBCn t•BC«i2l/X. 

C CtlCyiNU HC4T CONOUCtCD YNMUBN 01)H4CCB*Bf U/SIC 
I B4«B: I I BrN04*TBCN»4 l/X . 

IMCBNO:! .00)0 V*I64«B*4.T«) 1/ I )400.4l> . I 

ir IU4«B.0I.I040.I IHCO«0:i.))«4T«Y44«B*).«)XI/()4)0.*l/.l 
BCOMBO:4CONOe* INCOMO* I IBC 104- IBCNP* I /OCOOOO 
C CtUOLNIC MC41 CONOUCICO INNOUON riYCONAk. CM.B1U/)CC 
aCNNro.O 
DO )00 l:|*X 

IC4«4:|IC4Ntl I•IC4NI1•1II/X. 
tNC0N0:<.0O)BT4fC4VB*4.T9)l/«)4B9.*U.I 

IP IYC««t.BT.t040.l fNC0M0:l.0B«ftP*YC««B*S.9)7l/l>4D0.«IX.I 
)00 aC4N=«CM*4C*NI|l«YMCONONIICMIII-fC4Nlt*lll/IOC*NI)l*2.l 
C C41CUL4IC HCkl CONOUCUON INBOUBN BCBCHCNNIOO.BIU/Sf C 
tR4«6|:«)ttO*IRtl/X. 
tlia«SX:(lt|*IR2l/X. 

INCNOI = I.O))BI*IB4VBI*4.T9)l/l)4)ll.«tX.I 

IP I 194901.01. 1040. I fHCNOI:<.00*4l«l94V6l«$.9)n/l)400.*IB.I 
TNCNOXSI. 00 )aTN|R490I*4. 19)1/1)4)0 .NIX. I 

IP nB490X.oi.i040.i iNCNOx:i.oo«4i*iN4yex.s.9)xi/i)6ao.*ix.i 
0CRlNr).|»X*«RO**X-R|)9*2l*IMCN0l •IIRO-IRII/MLI 
8CROUl:4.Xa)»lNCNOX9RI)NIR)-ROI*IYRO-IN|l«IIICX*4tOOIMR|-RIII* 
IIRO*RI1II/IIR|*RI1I*IRO-R||IIII 
)CRio:iar«iN*«cRouYi/x.*a. 

C C4LCU14IC SMUTUC MC4Y IR4NSPCR. BIU/SCC 

tHCN00:2.«ait-O9N«TC*llll*ICni)ll/l.*l...; -04 

IF <NM04l.C0.NI tHCN00:|.94XC-09NI ICVK 1 1 • f C«U )l I/?.* I . I t»t -4 
OSMm:fHCNOON).INt4*01)PO*SIRORCN*?NIICnill-IC«l<)M/<a.O*CU4RN 
im)*iLN II 

aCIURN 
CNO 


Subroutine XOELP 


SUaaauMNt K0CIF<N.4RC41N,4RC40I .ONa.IlGtMtROl.VlOl.PlOMlN.PlINt 
ICOtP.OP.NI 

c IOC9IIPIC4IION or irn or los&rs 

C MTPin rONTR4CtION--*flPt:X CIP4NSI0N--MVPt:) 90-IUNM — KtVRCN I 

e Nl9*|:% OIHCR lUBN --UMPtS* PR Id t ON- TURI --NT VPt : I P R Id lON-RP bf N( 

C RtVRt:a PRICTtON*NONlMuN FOR lUBC — I SOiMt «N«L 
C m»RC = 9 PRICIIONMONCNIun POR RCOCNCR4fOR--lSOlNtRM«t 
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oi<ic«sto« cocriu%o>»OPK<«oi 

0«T« ••«*/]<.>/• I*SI/O.OMV«»/ 

•MNkSO* 

tr tm«>i»MtN».i.c.o«o> 0C1UBH 

p»o. 

B«uf:o. 

»P<lOilsO. 

aivBCsN 

•IdSABCalM/IBB. 

ttCtpHS/lt* 

Biscsvi%a*tir* 

riON:mtFio»tiii#a.i 

•OBCMSB* 

iFi»iit>»»wn sotMdi 
so kMINStlll 
•stlO/kOUT 

so to St 
SI kOlltsOBOl 
•slOOt/klN 
S? COOTIIIVC 

IFtk.lt. O.OOIt ksO.O 
• ClSFLOH/kHIII/llO 
»ClN0:00kVei««2/ <2.0kCB«« I 
BCSOO«*Ck*Ot/«lSC 

SO to lt•s*s.••».o.t•ott•»•Tlfktvn 
c 

1 CPOtlNUt 

C ktVPtSI COOtOkCtlON 

IFtttC-SOOO.QI 
n COCF:-3.k*k*3.S 
SO to IDO 

It COCFsHB.k*k«I.O 
SO to too 
c 

t COOtlNUC 

C NtVFC:? riFkMSION 

iMBC'SOOo.oi n,ti,tt 

tt COtF:il.9-?.k*k«l.}3S*k**SI*l-1.9l 

so to 100 

tt C0CF:ll.0>F.0«00*k*0.«90*««*2l*l*’l.0l 

to to too 

c 

s coktiNur 

CoeFsO.22 

so to tot 

c 

« CONtIMUF 

C0tF:S.0*0.2? 

so to 10 ? 

c 

s coktikur 

COtt £0.11 

60 to 10 ? 

c 

b CONItNUt 

C HttFfsb r«tCTI0N>1USt 

iriOC.lt.lSOO.i COCF.ib.O/OCklllBFlOt/k.OM 
IFIOE.6C.IS09.I C0CF:S.9kbF«t*O.?*tUS/l0e/«.0ll 
60 TO 10? 

C 

7 COkTlkur 

C RttPtSt rktCTI0ll>kE6EI«t»*t00 
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coer:o*«MC>ipi*o.oi«o*aci*o.M 
cotr:cocr«(iiis/iot/**oii 
to 10 102 

1)0 eotr:)*)*«**<*cocr 
10 ) DO:(OCr««rLNO*Otl 

in M«oc%LT«ai so io loi 

c <rvoc:o r*teiioii*No»it<tioN >>1001 «iv»c:« s*Ni 100 orstNfosie 

c mvotsi) ••1WSC NON.*r«u. kivk:ii*<-fo» otsr**ioo toi 
tfUTVOC.OC.lOl SSHHSISSSMI* 
IHMsvri/SOOnOOMNOOOlNolOS.O/OOOMkOI 
IH»):I.D/$0f>ni»O/NNt|**)-BkNH»«COiri 

rHS)i:i.o 

irussmH»t>iiH«ii/RN«ii.Lr.).})n ip 10 %t 

}| lN*0:llN«i*lNk)l*0.t 

rMk}|:MS«NN»-l.ai«>H*3**)*).)l/l tS*N»'*>IOI*RN»l««)*).OI 
t»lNlV0t.L1.t0l FHkfUl.O 
iriRtVPt.LldOt SkNNklSt.O 

knk):|.0)S001 lt•0/«Nk)*•)-Skl)Nk•COCr*Bk•l^k•lCkMHkl*l•Ol/).o/ekM••kl 
l•kkOGIIIIOkmNkOI••)•FHk^t II 

tnkOSIIitkf-RMkOl.Sl.O.OOl I so '4 SI 
S) m:i.)/Fiik)l 

«o:xi«k)/kNkt«soiitii.o/rNk)ii 
»IISINkt/l|Nk2«$Q01l|.0/rMk)|l 
»0(l):PllkPlN 
OP:Pltt*»OUT 
1)1 COSIINUt 

inPkOMIO.lC.O.OI OP:|-l*l«OP 
COCPRINISCOCP 
OPXIkIsOP 
ktiuilii 
CNO 
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Printed Output ol Sample Run 


aamaa 

e 

s 

.aagBSBBac*9S 




auaat 

r 

•tSSSS9B3t*>l 




riecac 

t 

•IBOBBBBB€«S« 




teaiat 

t 

•SSB 




ttnea 

s 

•a 




cacPT 

s 

• IBOB9BOBC-BS 




IS 

: 

»Togg99ggf*BS 




es 

s 

•anssB03t*ss 




oacaa 

X 

•aOBBaOBBtOS 




racTi 

: 

•aa8BB0B8t*BB 




raeiT 

X 

•ISSSSS6SC*3f 




aacic 

s 

•SB 




aitai 

s 

•a 




aaeao 

X 

•as 




Pages 

3 

*a 




ta 

3 

..iao9oaBC*aa* 

•laaaoooetPBa. 

•iaa90B9BC«8a* 

.iaaDoaooc«j%i 



.taaas9osr«sa. 

.ina«0B09C*aa* 

•tS9a999BC«9a« 

.ior«Ba9B(«Qai 

ran 


•TT|PSpg9t*89« 

•a%6goooof*os. 

-eamigaaot^os. 

.a«BOBBf<OC*OI 

S 

• laBBBBOBCOa* 

•iaaoooooc*o«. 

• iaaoBeao(*Ba • 

•iaeoBoaoc«oa, 



•lasaaoBMps** 

•iaa«B8B8t»9a* 

•tS9aS8SBt*Sa* 

• ia?«B399C*0%( 



• TTSOSBBBI^BS* 
•aasssaoBc*ss 

•%aaooooBC*os« 

• a%0P80S0t*B9. 

•aaoBeeooc«09, 

in 

t 

• lTaBSBOBC«aa. 

.iasooBOOt*o%. 

•iaS80O9OC*Ba. 

•i6Ssoooor*oa( 



•iasassBSt*9«. 

.iaa%B909C*9a* 

•tSSaS399€*9%. 

•l9T«B390C*Oa, 



•TTa9S009C*BS. 

•aaoppooBtos 

•saQoooogr*gs. 

•SaOBBOBBC^DS. 

•aaBOBgOBC«P», 

anCMC 

S 

.IOB09000t«St 




Jt» 

~ 

•I 





acefi 

laat c»ut 


ales sio.tM aaiiTs i?e.»u eentkpr soa.elT eancaes-tai.}** amiTof: freietr «%m acrri: •eer acrf?: 

o»C6Ca: . 1 ST eeaep: is.isr ferg: aa.sss uMan iT.ata chims -i.ass uieiH: ifgo.ees 

usee: ist6.aa« avaase: tto.ud gaiae: -iB.sit iieisi: laa.aoo ateisi: isa.noo «ene: a.aas 

Miens -as.wa Uaees at.ara gNraies-ssa.eiK one ate: so.aa? ocsoias .bop ocoeir: is».Taa ecoeea: -.pst ecoeee 

OPCSn: -asum BBfaes ast.ieo HDisees a.TiB aot&efis -ia.au «eisies sai.aia eeiains-jaT.tu oupeis a.aas acnoeo 

BCapCis a*taT aendOs i.sss ocasens i.soi TacieasiTTa.iTa Tacnea: asa.ooa qsnTi: i.aaT 

oi«a s sia.76t HPaaaa s lai.aaa 

Miaaat mats a» tasT t enus 


ains ssa.aa? aoMts ira.m eaarie: i9a.a«a aaacees-iap.tTa iimTai: iis.sia crMprs .«S6 errrit 
oecau: .us eaeae: is.sfa fecps aa.DM aeaiis if.aai ruins -i.aai utouls irta.aaa 

usee: isaa.aaa aaaase: TaT.ear nptaes -lo.us neisis laa.ps? uTeist: tsa.Dri nmiLe: a.aj? 

Buens -ss.aaa acaee: «).ari enratas-ssa.saa one ties la.aas ocaoLas .jsp otooies laa.TQv oconpar 
gacans *a?a.«os aacae: ara.sai noispe: a.aaa aoiseas •la.aaa ppisips aia.avo aemns-ias.aaa gmoais 

acBBCLs a.iar acnoos i.asa oenoens i.ssi lacieanm.STa tacneas an.aao ganTi: «,«aT 

aina : m.ias uaaaas = lai.sgg 
aataaar ttnetaaivati o«ca list ctcic 

vaoes iTOa* iti«. lass. laaa. tsan isoa. lois. aii. eat. es«. eta. 

laacTc: m«. it6p. i?*}« lara. i«ar. ms. tiia. 9?«. t* 5 . esi. au. tu 

mcTCs lua. lasg. iiso. u;a. i«as» ups. itis. tti. ns. &«o. seo. s«o 


aai airr?s 


-.OSS QCOnep 
i.aat octtoeo 


a*«. 

• ate. 

. s«o. 


.♦9P 

2.SOP 

: a*«ss 

T.Tia 

I.pit 
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DEFINITIONS OF NAMELIST INPUT VARIABLES* 


COEFF leakage coefficient for flow between buffer and 

compression spaces 

FACTl, FACT2 ccunvergence constants used in regenerator metal 

temperature convergence scheme 

FIPCYC number of Iterations per cycle 

IPRINT number of iterations between printouts (if JIP=0) 

TTMPS 1, printout gas and metal control volume temper' 

ature at every iteration 
0, do not do the above 

JIP >0, short •form printout 

=0, long-form printout 

MWGAS 4, use helium working gas 

use hydrogen working gas 

NOCYC number of engine cycles to be run 

NOEND number of cycle at which convergence scheme is 

turned off 


NSTRT number of cycle at which regenerator metal temper- 

ature convergence scheme is turned on 

C^EGA drive speed, Hz 

2 2 

P initial pressure, Ibf in, (N/m ) (approximately 

equal to mean pressure for GPU) 

P3 initial buffer space pressure. Ibf /in. (N/m ) 

REALGS 1, use real -gas equation of state 

0, use ideal -gas equation of state 


^Although both U.S. customary and SI units ai'e given here, U.S. 
('ustomary units are used in the prograin. 
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RHCFAC regenerator heat -transfer coefficient multiplication factor 
TG 12 control volume interface gas temperatures, ^ ( K) 

TGA 13 control volume gas temperatures, ^ ( K) 

TM 13 metal boundary temperatures, ^ ( K) 

T3 buffer space temperidure (assumed constant), °1R ( K) 

D^nitions of Output Variables 

AVGWSP mean working -space pressure, Ibf/in. (N/m ) 

EFFTOT efficiency, WRKTOT/QOJ 

ENTH net enthalpy flow from working space to buffer space over 
one cycle, ft-lbf (J) 

O 

EN3PU rate of displacement of working -space gas, in. /min 
(cm^/mta) 

FREQ relational frequency of craidc^aft, He 

PWRHP Indicated power - WRKTOT/cycle period, hp (kW) 

QCNIX7L heat conduction through cylinder, ft-lbf /cycle (J/cycle) 

QCNDCN heat conduction through insulation container, ft-lbf/cycle 
(J/cycle) 

QCNDD heat conduction through displacer, ft-lbf/cycle (J/cycle) 

QCNDRI heat conduction through regenerator casing, ft-lbf/cycle 
(J/cycle) 

QCNDRO heat conduction through r^enerator casing, ft-lbf/cycle 
(J/cycle) 

QCOMPN heat from compression-space wall to gas, ft-lbf/cycle 
(J/cycle) 

QCOMPP heat from gas to compression -space wall, ft-lbf/cycle 
(J/cycle) 
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QCOOLN heat firom c<M>ler to gas, ft-lbf/cycle (J/cyole) 

QCCX>LP heat from gas to cooler, ft -Ibf/ cycle (J/cycle) 

QEXPN heat from metal wall to gas in expansion space, ft*lbf/cycle 
(»T/cycle) 

QEXPP heat from gas to metal wall in expansion space, ft-lbf/cycle 
(J/cycle) 

^EATN heat from heater to gas, ft-lbf/cycle (J/cycle) 

QHEATP heat fbom gas to heater, ft-lbf/cycle (J/cycle) 

QIN heat into gas in expansion space ami heater tubes over one 

cycle, ft -Ibf (J) 

QOUT heat out of gas in cooler tubes and compression space over 
one cycle, ft -Ibf (J) 

^EGEN net heat flow from gas to regenerator metal over one cycle, 
ft -Ibf (J) 

QREGN heat brom regenerator metal to gas, ft-lbf/cycle (J/cycle) 

QREGP heat from gas to regenerator metal, ft-lbf/cycle (J/cycle) 

QSHTL shuttle heat loss, ft-lbf/cycle (J/cycle) 

REFFl a measure of reg.''nerator effectiveness 

REFF2 another measiure of regenerator effectiveness 

TG gas temperature at interface of two adjacent control volumes, 

®R ( K) 

TGA working bulk gas temperature of a control volmne, ^R ( K) 

TGACYC cycle-time average temperature of TGA, °R ( K) 

TGCMPA as TGACYC in compression space 

TGCYC cycle-time average temperahme of TG, °R ( K) 


TGEXPA as TGCYC in expansion space 

TM metal temperature of a metal control volume, ®R ( K) 
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TMCYC 

UTOTAL 

WLiSP 

WDBPN 

0WDISPP 

WPISTP 

WPIST 

WPISTN 

WRKCMP 


WRKEXP 

WRiJp 


WRKLT 

WRKl’OT 


cycle-time a\ erage temperature of TM, ( K) 

internal energy content of working -space gas, ft-lbf (J) 

net work done on displacer piston by working -space gas, 
ft-lbf (J) 

work done on gas by displacer, ft-lbf /cycle (J/cycle) 

work done on displacer by gas, ft-lbf /cycle (J/cycle) 

work done on piston by gas, ft-lbf /cycle (J/cycle) 

work done on power piston by workii^ -space gas, ft-lbf (J) 

work done on gas by piston, ft-lbf /cycle (J/cycle) 

work done by tlie ^as in compression space over one cycle, 
ft-lbf (J) 

work done by the gr^ in expansion space over one cycle, 
ft-lbf (J) 

power -piston work lost due to pressure drop from regen- 
erator to compression space, ft-lbf /cycle (J/cycle) 

^rk lost due to pressure drop over one cycle, ft-lbf (J) 

net work done by gas in working space over one^^cle, 
WRKEXP + WRKCMP, ft-lbf (J) 
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TABLE I. - CHARACTERISTICS OF STIRLING ENGINE 


GROUND-POWER -UNIT SIMULATION RUN 


Working fluid Hydrogen | 

Frequency, He 50 ■ 

Average vrorklng -space pressure, N/m^ Opsl). . . 5. 50x10® (800) j 


Heater temperature (average gas temperature for i 

coidrol volume 2, fig. 2), K (®R) 976 (1760) | 

Cooler temperature (cooler metal -wall temperature), 

K (^) 300 (540) 

Leakage flow coefficient (used in eq. (27)) 0. 0001 


TABLE n. - PERFORMANCE PREDICTIONS FOR STIRLING ENGINE 
GROUND-POWER -UNIT SIMULATION RUN 


Indicated power, kW (lq>) 11. 0 (14. 8) 

Indicated efficiency 0.49 

Indicated work per cycle, J (ft-lbf) 220.4 (162.5) 

Work loss per cycle, due to pressure drop, J (ft-lbf) 1 1. 7 (8. 6) 

Conduction heat loss per cycle, J (ft-lbf) 20.9 (15.4) 


j Shuttle heat loss per cycle, J (ft-lbf) 8. 3 (6. 1) 




Variable 

Metai vconstantU temperaturj^ Constanj^ 




Inputs: Frequency of qieration 

Pressure level (of initial voiumel 
Initial temperatures 


Revise 

regenerator metal 
temperatures 
to speed up 
convergence 


No I Yes 
18. Revise 
regenerator 
metal 

temperatures? 


No 


L 

2. 

3L 

4 

% 

4 

7. 


la 

IL 

IZ 

131 

14 

15. 


=n 

Tine - Time 4(Time) 

Update pressure Imml. P - f (volumes, temimratures. massesi 

Update teiiqieratur*s In each gas control vtduim for effect of change 
in ^assure levd 

Update nuss distriiNJtion 
Update flow rates 

Uptete temperatures in each gas cmitrol v(dume for effect of flow 
between control vfdumes 

Upd^ heat-transfer coefficients between motel and gas control 
volumes 

Update temperature in each gas contrd volume for effect of heat 
transfer (between gas and meted 

Update regenerator metal temperature for effect of heat transfer 
(between ge and meted 

(telculate presure drops (neesary over only one cycle) 

Calculate conduction and shuttle lese (just once during last cycle) 
Sum up heat transfers between ge and metal far ech component 
Uptete iekage flow b^ween working and buffer space 
Update expansion-, cmnpresimi-. and buffer-space volume 

Calculate work done in expansion and compresion space and sum 
up to get totei indicated work 


No 


15. is cyde complete? 

I Ye 

14 Calculate indicated iXMer and dficiency 
17. Have specified number of cycie been completed? 
Ye 

Call Calcomp to plot engine variable if requeted 


Figure 3L - Outline of caiculation procedure. 
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a aova Tvinow^’ 





Temperature 



(a) Sample regenerator control volume: 



I-l I 

Regenerator node 


(b) Control-volume temperature profile. 

Figure 7. - Sample regenerator control volume and temperature profile. 



Points derived 
from data 



Figure 8. - Regenerator friction factor data and curve fit 
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Figure 9. - Schematic showing dimensions needed for calculating heat conduction. (Regenerator, housing, cylinder, 
and displacer are 310 stainless steel. Dimensions are in cm (in. I. ) 





Figure 10. - Schematic showing geometric relations between piston 
positions and crankshaft angle. 




AAetei boundary 
temperature 



Figure IL - Assumed metal boundary temperat^^' < for sample run. 
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Time, sec 

Figure IZ - Expansion' and compression -space wiumes as functions of time. 
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Time, sec 

Figure 15. - Pressure drop as a function of time. 



Figure 16. - Heater flow rates as functions of time. 
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Compression-siiuca end 



Dd$/qi m\i jaiooQ 


Figure 18. - Cooler flow rates as functions of time. 



0 1 2 3 4 5 6 7 8 

Expansion-space tfolu me. in^ 

r v''e 1% - Expansion-space pressure-volume diagram. 


Compression-space pressure, Ib/in 


8. Si 


Csl 




Figure 2d - Compression-space pressure-voiume diagram. 


Pressure, 



0 2 4 6 8 10 12 14 16 18 20 

Working-space volume, in^ 


Figure 2L - Total-wortung-space pressure-volume diagram Ineglects pressure drop). 




§ ^ P n 


X *ajn)ej9dui9) 99eds-uo!ss9Jdujoo 



Hq *9jn)9J9duj9) 99eds-uo|ss9Jduioo 
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1 

2 

20 30 40 50 60 

Engine frequency, Hz 

Figure 24 ' Indicated power as a function of engine frequency for 
several average working-space pressures. 






Indicated efficiency, percent 



Figure ■ Indicated efficiency as a function of engine 
freni'f -cylpr several average working-space pressures. 
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Engine frequency, Hz 

Figure 27. - Indicated efficiency as a function of 
and helium. 


Engine frequency, 



)Ud3jdd 'Aouepiue pe^eoipui 


Figure 28. - Indicated efficiency as a function of indicated power (with engine 
frequency shown on curves). 



Figure 29. - Comparison of equations (C3) and (C4). 



N 'Jaqiunu Jiassn|\| 


Figure 30. - Regenerator heat-transfer data used for ground- 
power-unit model. 
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